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ABSTRACT 


The  progressive,  sinusoidal  water- wave-induced  perturbations 
in  a turbulent  air  flow  above  an  air-water  interface  were  studied. 
The  velocity  field  was  measured  both  in  fixed  and  wave-following 
reference  frames  with  a cross  hot-film  probe.  Both  the  magnitude 
and  phase  of  the  wave-induced  velocities  and  turbulent  Reynolds 
stresses  were  obtained  across  the  turbulent  boundary  layer. 

The  mean  momentum  transfer  from  the  wind  to  water  waves  is 
comprised  of: 

(1)  A critical  layer  part,  which  is  proportional  to  the 
product  of  the  velocity-profile  curvature  and  the  mean  square  of 
the  wave-induced  vertical  velocity  at  the  critical  height,  where 
the  mean  wind  speed  is  equal  to  the  wave  speed. 

(2)  A vertical  integral  across  the  turbulent  boundary  layer 
of  the  mean  product  of  the  wave-induced  vertical  velocity  and 
vorticity  associated  with  the  wave-induced  turbulent  Reynolds 
stresses . 

(3)  A contribution  associated  with  the  wave-induced~tur- 
bulent  Reynolds  stresses  at  the  air-water  interface. 

These  quantities  were  evaluated  for  different  wind  speeds. 
The  calculated  results  showed  that  the  wave-induced  turbulent 
Reynolds  stresses  are  important  in  the  evaluation  of  momentum 
transfer  from  wind  to  waves,  their  contribution  comprising  about 
90  percent  of  the  total  transfer,  while  the  critical  layer  con- 
tribution is  about  10  percent. 
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1.  INTRODUCTION 


1 . 1 Motivation  and  Previous  Work 

The  search  for  an  understanding  of  the  physical  mechanisms  in- 
volved in  the  generation  and  maintenance  of  waves  by  wind  has  been 
one  of  the  most  difficult  in  the  field  of  fluid  mechanics.  It 
is  only  in  the  past  two  decades  that  much  progress  has  been  made. 

In  a survey  paper,  entitled  "Wave  Generation  by  Wind",  Ursell (1956) 
commenced:  "Wind  blowing  over  a water  surface  generates  waves  in 

the  water  by  physical  processes  which  cannot  be  regarded  as  known". 
His  review  stimulated  wide  interest;  since  then  many  investigators 
have  contributed  to  the  clarification  of  the  processes,  and  a 
general  theoretical  frame-work  has  emerged.  The  first  significant 
contributions  to  the  dynamical  wave  theory  were  Phillips'  (1957) 
and  Miles'  (1957)  theories,  which  yielded  rigorous  transfer  expres- 
sions for  certain  aspects  of  the  wave-atmosphere  interactions. 

Phillips'  (1957)  resonance  mechanism  accounts  for  the  initial 
linear  growth  of  a wave  spectrum  under  the  influence  of  random  pres- 
sure fluctuations.  The  statistical  properties  of  the  pressure 
fluctuations  are  associated  with  the  atmospheric  turbulent  wind 
field  and  assumed  to  be  independent  of  the  waves  generated.  The 
assumptions  fail  when  wave  amplitude  is  sufficiently  large.  If 
the  waves,  having  already  been  initiated,  are  of  sufficient  ampli- 
tude the  air  flow  over  the  undulatory  surface  will  cause  the 
generation  of  pressure  differences  along  the  perturbed  water 


surface.  Miles'  (1957)  inviscid  model  described  the  energy  trans- 


f 


fer  between  the  wind  and  the  water  wave  by  the  wave-induced  pres- 
sure component  that  is  in  phase  with  the  vertical  velocity  com- 
ponent of  the  water  surface.  The  model  considered  air  flow  over 
a single  frequency  progressive  wave  of  infinitesimal  amplitude. 
Turbulence  in  the  air  stream  was  neglected  except  in  maintaining 
a prescribed  parallel  shear  flow  with  logarithmic  velocity  distri- 
bution. The  feedbacks  from  the  wave  were  the  wave-induced  air 
velocity  and  the  pressure  perturbation.  The  wave-induced  air 
velocity  was  assumed  to  be  two  dimensional  and  small  enough  to  be 
unimportant  in  the  nonlinear  processes  of  the  equations  of  motion. 

The  flowfield  in  the  water  was  considered  inviscid,  irrotational 
and  incompressible.  Mean  water  currents  were  assumed  to  be  absent. 
The  equations  governing  the  wave-induced  perturbations  in  the  air 
stream  are  exactly  the  same  as  those  in  laminar-instability  theory. 

Miles'  inviscid  model  has  been  tested  both  in  the  field 
and  in  the  laboratory  through  the  measurements  of  wave  growth-rate 
and  wave-induced  pressure  in  the  air  stream.  Wave  growth  measure- 
ments in  the  field  by  Snyder  and  Cox  (1966)  and  by  Barnett  and 
Wilkerson  (1967)  indicated  that  the  observed  exponential  growth 
rate  exceeds  that  predicted  by  the  inviscid  model  by  an  order  of 
magnitude.  The  wave  growth-rate  obtained  from  Dobson's  (1971) 
power  and  cross-power  spectra  of  wave  elevation  and  pressure  agrees 
with  the  field  measured  wave  growth-rate,  but  the  measured  wave 
growth-rate  is  larger  by  a factor  between  5 to  8 than  Miles'  inviscid 
theoretical  predictions. 

Shemdin  and  Hsu  (1967)  investigated  the  wave-induced  pres- 
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sure  distribution  over  mechanically  generated  waves.  Their  re- 


sults indicated  that  the  growth  rate  was  larger  by  a factor  of 
about  two  than  that  indicated  by  the  inviscid  results.  Bole  and 
Hsu  (1969)  studied  the  growth  of  mechanically-generated  waves 
under  the  action  of  wind  in  a channel.  Based  on  the  measured  spa- 
tial growth  of  waves  along  the  channel,  they  obtained  the  growth 
rate  of  the  waves  and  found  the  measured  growth-rate  larger  by  a 
factor  from  one  to  ten  than  Miles'  estimates.  In  order  to  examine 
the  effect  of  neglecting  the  viscosity  in  Miles'  inviscid  model 
(the  quasilaminar  model),  Benjamin  (1959)  and  Miles  (1959)  in- 
cluded viscous  term  in  the  equations  of  motion  and  carried  out  the 
analysis  in  curvilinear,  orthogonal  coordinates.  They  found  that 
the  viscous  effects  are  of  second  order  and  can  be  neglected. 

A further  contribution  to  dynamical  wave  theory  was  made  by 
Phillips  (1960)  and  Hasselmann  (1962,  1963).  They  studied  the 
energy  transfer  due  to  nonlinear  wave-wave  interactions.  The 
most  important  simple  result  of  these  calculations  is  that  the  time 
scale  for  the  interaction  is  much  greater  than  the  time  scale  of 
wave  growth;  the  wave-wave  interactions  cannot  transfer  energy  from 
a band  of  wave  number  as  rapidly  as  it  is  acquired  from  the  wind. 
Moreover,  the  energy  transfer  among  the  different  components  is 
weak  and  selective.  Only  for  certain  combinations  of  wave  numbers 
and  frequencies,  which  satisfy  the  resonance  condition  does  the 
nonlinear  interaction  produce  a continuing  energy  flux  from  one 
to  another. 
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Hasselmann  (1968)  applied  the  theory  of  weak  interaction  In 


1 

random  fields  to  describe  the  wave  field  within  the  coupled  ocean- 
atmosphere  system.  In  his  description,  the  lowest-order  of  trans- 
fer expressions  for  wave-atmosphere  interactions  are  the  Phillips  il 

and  Miles  processes,  a nonlinear  correction  to  Miles  process,  and 
wave-turbulence  interactions.  He  stated  that:  "Present  data 

suggests  that  the  wave-turbulence  interactions  may  be  the  most 
important  of  the  four". 

The  influence  of  the  background  turbulent  field  on  the  wave- 
induced  perturbation  appears  in  the  governing  equations  for  the 
wave-induced  field  as  an  oscillation  of  the  background  turbulent 
Reynolds  stresses.  The  significance  of  the  wave-induced  turbulent 
Reynolds  stresses  were  revealed  from  the  failure  of  Miles'  invis- 
cid  models  (1957,  1959),  in  which  the  wave-induced  turbulent 
Reynolds  stresses  were  neglected. 

In  predicting  the  wave-induced  velocities  and  pressure  in  the 
turbulent  air  stream,  the  difficulties  are  that  the  dynamic  equa- 
tions describing  those  quantities  in  a turbulent  flow  are  not  closed, 
due  to  the  presence  of  induced  turbulent  Reynolds  stresses.  There 
are  more  unknowns  than  equations.  Efforts  to  seek  closure  by  de- 
riving dynamic  equations  for  the  induced  turbulent  Reynolds  stresses 
cannot  solve  the  closure  problem  but  only  create  more  unknown  quan- 
tities. But  the  equations  for  the  induced  turbulent  Reynolds 
stresses  (see  Reynolds  and  Hussain  1972)  indicate  that  the  induced 
turbulent  Reynolds  stresses  are  of  the  same  order  of  magnitude 
as  the  other  disturbance  quantities  in  the  air  stream.  Thus,  even 
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for  weak  organized  motions  in  the  air  stream  one  should  expect 
oscillations  of  the  turbulent  Reynolds  stresses  of  comparable 
magnitude . 

Because  of  the  lack  of  an  adequate  theory  of  turbulence  at  the 
present,  plausible  assumptions  have  been  made  through  the  use  of 
the  turbulent  eddy  viscosity  and  turbulent  kinetic  energy  models 
in  which  induced  turbulent  Reynolds  stresses  are  retained.  The 
eddy  viscosity  and  turbulent  energy  closure  models  have  been  car- 
ried out  by  Davis  (1970),  Hussain  and  Reynolds  (1970),  Saeger 
and  Reynolds  (1971),  Long  (1971),  Townsend  (1972),  Davis  (1972), 
Acharya  and  Reynolds  (1975)  and  Norris  and  Reynolds  (1975).  Their 
results  were  compared  with  the  experimental  data  of  Stewart  (1970), 
Hussain  and  Reynolds  (1970),  Saeger  and  Reynolds  (1971),  Dobson 
(1971),  Acharya  and  Reynolds  (1975),  and  Norris  and  Reynolds  (1975). 

Davis  (1970)  proposed  that  the  Reynolds  shear  stress  of  the 
background  turbulence  is  distorted  by  an  amount  proportional  to 
the  displacement  that  streamlines  would  have  in  an  inviscid  uni- 


form velocity  flow  over  the  wavy  boundary,  the  proportionality  con- 
stant  being  the  normal  gradient  of  the  turbulent  shear  stress  in 


the  undisturbed  flow.  He  also  proposed  that  the  ratio  between  the 
different  components  of  induced  turbulent  Reynolds  stresses  is 
the  same  as  that  of  the  turbulent  Reynolds  stresses  in  the  turbu- 
lent boundary  layer  above  a flat  rigid  surface.  Numerical  pre- 
dictions based  on  this  model  are  inconsistent  with  Stewart's  (1970) 
wave-induced  velocity  data.  Generally,  at  high  wind  speeds  the 
predicted  magnitudes  of  wave-induced  velocities  are  larger  than  the 
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measured  values,  whereas  at  low  wind  speeds,  the  inverse  is  often 
the  case.  Davis  also  did  quasi laminar  calculations,  whose  pre- 
dicted values  of  wave- induced  velocities  showed  little  agreement 
with  the  experimental  data  of  Stewart  (1970);  the  amplitudes  tend 
to  be  too  high  at  high  wind  speeds  and  the  phases  show  poor  agree- 
ment. The  calculations  also  show  the  predicted  wave-induced  velo- 
city data  are  sensitive  to  details  of  the  mean  velocity  profile 
near  the  wave  surface  where  no  accurate,  reliable  experimental 
data  exist. 

Hussain  and  Reynolds  (1970)  studied  the  behavior  of  a periodic 
disturbance  in  turbulent  shear  flow.  They  postulated  that  the  in- 
duced turbulent  Reynolds  stresses  are  proportional  to  wave-induced 
strain  rate  through  an  eddy  viscosity  and  that  the  eddy  viscosity 
for  the  perturbed  flow  is  the  same  as  that  for  the  undisturbed  flow. 
The  eddy  viscosity  expression  across  the  turbulent  boundary  layer 
used  in  their  computation  was  the  one  proposed  by  Reynolds  and 
Tiederman  (1967).  The  relative  success  of  their  eddy  viscosity 
model  over  the  quasilaminar  model  in  predicting  wave  behavior  demon- 
strates that  the  distortion  of  the  background  turbulent  field  by 
the  wave  is  important  and  should  not  be  ignored.  Saeger  and  Reynolds 
(1971)  studied  the  interaction  of  sinusoidal  travelling  waves  with 
turbulent  shear  flow  in  a two  dimensional  channel  flow  with  one  me- 
chanically articulated,  waving  wall.  The  turbulent  closure  assump- 
tions for  the  wave  perturbed  momentum  equations  are  same  as  those 
of  Hussain  and  Reynolds,  namely,  the  eddy  viscosity  model.  The 
measured  oscillating  pressures  were  compared  with  predictions  of 
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inviscid,  quasilaminar  and  turbulent  eddy  viscosity  models.  They 
found  the  trends  of  pressure  amplitude  and  phases  are  correctly 
predicted  for  upstream  running  waves  by  all  three  models.  For 
downstream  running  waves  the  models  begin  to  differ  substantially, 
namely,  neither  the  amplitude  nor  phase  agree  with  the  data.  Also, 
the  inviscid  model  predicts  smaller  values  of  power  transfer  than 
those  predicted  by  the  turbulent  eddy  viscosity  model. 

The  eddy  viscosity  model  implies  that  the  effect  of  turbulent 
mixing  is  similar  to  the  effect  of  molecular  viscosity  in  that  the 
stress  is  proportional  to  the  velocity  shear.  The  model  also  im- 
plies that  the  turbulent  Reynolds  stresses  is  a local  phenomenon; 
its  effect  at  each  point  is  determined  by  the  character  of  the  mean 
flow  at  that  point.  Noting  that  the  flux  of  momentum  which  the 
Reynolds  stress  represents  is  largely  the  result  of  mixing  between 
regions  where  the  mean  flow  differs,  Prandtl  (1945)  suggested 
that  the  turbulent  Reynolds  stresses  must  be  related  to  the  energy 
of  turbulent  fluctuation.  Long  (1971)  adapted  the  turbulent  energy 
closure  model  of  Bradshaw,  Ferriss  and  Atwell  (1967)  to  obtain 
similar  closure  relations  for  wave-induced  fields.  The  results  of 
his  model  were  compared  with  the  results  of  Miles  inviscid  theory 
and  the  experimental  work  of  Dobson.  The  experimental  results  of 
Dobson  agree  better  with  predictions  of  Long's  (1971)  turbulent 
model  than  with  Miles'  inviscid  theory.  In  some  typical  cases  the 
computed  growth  rate  differs  from  Dobson's  experimental  values  by 
a factor  no  more  than  two.  On  the  whole,  the  comparisons  between 
experimental  and  theoretical  growth  rate  spectra  and  pressure 
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phase  angles  showed  better  agreement  at  the  lower  wind  speeds  than 
at  the  higher  wind  speeds.  It  was  also  reported  that  the  computa- 
tion was  very  sensitive  to  the  boundary  condition  of  wave-induced 
turbulent  energy  at  the  air  water  interface. 

The  visco-elastic  behavior  of  a turbulent  fluid  was  added  to 
the  Bradshaw,  Ferriss  and  Atwell  (1967)  scheme  by  Townsend  (1972) 
to  investigate  the  distortion  of  turbulent  boundary  layer  flow  due 
to  surface  waves  running  in  different  directions  from  the  mean 
velocity  field.  The  induced  turbulent  Reynolds  stresses  were  charac- 
terized by  elastic  elongation  and  rotation  of  the  mean  turbuLent 
Reynolds  stresses.  The  calculated  wave  growth  rates  are  consi- 
derably less  than  those  measured  by  Snyder  and  Cox  (1966),  Barnett 
and  Wilkerson  (1967),  and  Dobson  (1971). 

Both  the  turbulent  kinetic  energy  and  eddy  viscosity  models 
were  reconsidered  by  Davis  (1972).  He  used  a phenomenological 
closure  of  the  conservation  equation  for  turbulent  Reynolds  stresses 
to  test  two  generalizations  of  the  Bradshaw,  Ferriss  and  Atwell 
(1967)  model.  In  one  all  stresses  are  assumed  to  vary  proportionally 
to  mean  turbulent  energy;  in  the  second  only  the  normal  stress 
variations  are  proportional  to  turbulent  energy  and  the  conservation 
equation  of  the  induced  horizontal  shear  is  retained.  The  computed 
growth  rate  for  the  first  generalization  was  an  order  of  magnitude 
less  than  the  experimental  results  of  Kendall  (1970)  and  Dobson ( 1 971 ) . 
The  computed  results  for  the  second  generalization  were  worse. 

Davis’  work  also  made  clear  that  the  models  are  very  sensitive  to 
the  mean  velocity  profile  near  the  air-water  interface.  This 
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makes  it  difficult  to  draw  any  definite  conclusions  from  com- 
parisons with  experiments  since  the  mean  velocity  profile  is  un- 
known near  the  wave  surface. 

Norris  and  Reynolds  (1975)  studied  the  wave-induced  pressure 
in  a two  dimensional  channel  flow  with  one  mechanically-articulated 
waving  wall.  The  pressure  measurements  were  on  the  stationary  wall 
opposite  to  the  moving  wavy  wall.  They  solved  a mathematical  model 
of  the  problem  in  a transformed  coordinate  system  so  that  the 
boundary  conditions  are  applied  exactly  at  the  waving  wall.  They 
reported  that  their  turbulent  kinetic  energy,  eddy  viscosity,  and 
quasilaminar  models  produced  about  equally  good  pressure  predic- 
tions; the  use  of  the  coordinate  transformation  seemed  crucial  in 
their  study  and  the  choice  of  Reynolds  stress  model  is  not  impor- 
tant. They  found  distinctly  better  results  with  the  transformed 
coordinate  system  (recall  Benjamin  (1959)  and  Miles  (1959)  also 
used  similar  curvilinear  coordinates).  Acharya  and  Reynolds  (1975) 
used  the  same  turbulent  model  as  Norris  and  Reynolds  (1975)  to 
study  the  response  of  turbulent  channel  flow  to  imposed  oscillations. 
Their  major  conclusion  is  that  the  dynamics  of  the  turbulence  must 
be  considered  in  any  prediction  of  unsteady  turbulent  flows  and  a 
Reynolds  stress  equation  model  must  be  used. 

Besides  the  phenomenological  turbulent  models  mentioned  above, 
Davis  (1974),  working  on  the  basis  of  several  assumptions,  used 
a non-phenomenological  model  to  investigate  the  generation  of 
Reynolds  stresses  by  an  infinitesimal  perturbation  of  a turbulent 
shear  flow.  The  general  method  is, in  principle,  capable  of 
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producing  a constitutive  relation,  but  unfortunately  putting  his 
method  into  practice  requires  a detailed  description  of  the  undis- 
turbed turbulent  field  which  is  considerably  more  comprehensive 
than  anything  that  can  be  obtained  from  presently  available  experi- 
mental data. 

The  possible  significance  of  the  induced  turbulent  Reynolds 
stresses  on  the  momentum  transfer  from  wind  to  wave  was  pointed  out 
both  by  Miles  (1967)  and  Phillips  (1966).  Kendall  (1970)  experi- 
mentally studied  the  response  of  the  turbulent  flow  structure  to 
the  perturbation  imposed  by  a wavy  wall.  He  indicated  experimentally 
the  strong  modulation  of  turbulent  structure  due  to  progressive 
waves  and  the  non-negligible  energy  transfer  to  the  waves  due  to 
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induced  turbulent  Reynolds  stresses. 

Yu,  et  al.  (1973)  attempted  for  the  first  time  to  measure 
the  induced  turbulent  Reynolds  stresses  over  a progressive  water 
wave  in  the  turbulent  boundary  layer.  Because  of  the  limited 
accuracy  provided  by  analog  techniques  in  data  taking  and  analysis 
their  results  were  only  qualitative  and  unsatisfactory. 

1 . 2 Obj ectives  of  the  Study 

The  experimental  program  was  carried  out  in  the  Stanford  Wind, 
Water-Wave  Research  channel  with  the  objective  of  studying  the 
interactions  between  a turbulent  shear  flow  of  air  and  a progres- 
sive, sinusoidal  water  wave.  Both  the  magnitude  and  phase  of  wave- 
induced  velocities  and  turbulent  Reynolds  stresses  across  the  tur- 
bulent boundary  layer  were  obtained. 


The  specific  objectives  of  this  study  are  as  follows: 


a.  To  reveal  the  nature  of  both  the  magnitude  and  phase  dis- 
tributions  of  wave-induced  velocities  and  turbulent 
Reynolds  stresses  across  a turbulent  boundary  layer.  With 
a newly  acquired  digital  data-acquisition-reduction  sys- 
tem, we  hoped  to  obtain  experimental  data  of  greatly  im- 
proved quality  and  valuable  information  to  facilitate  the 
development  of  improved  closure  models. 

I 

b.  To  evaluate  the  importance  of  turbulence-wave  interactions 
in  the  formulation  of  wind-wave  generation  theory.  Spe- 

I cifically,  from  the  measured  wave-induced  velocities  and 

turbulent  Reynolds  stresses,  the  importance  of  turbulence 
effect  on  the  momentum  transfer  from  wind  to  waves  could 
be  verified. 

c.  To  clarify  the  different  and  contradicting  conclusions 
reached  by  various  investigators  regarding  the  importance 

| 

of  the  wave-induced  turbulent  Reynolds  stresses  in  rela- 
t ■ tionship  to  different  techniques  employed  in  their  analyses. 


ties.  The  relation  of  the  wave-induced  turbulent  Reynolds 
stresses  to  wind-wave  generation  theory  is  also  examined. 

b.  In  Chapter  3,  the  experimental  facilities  and  instrumen- 
tation used  to  study  the  interaction  between  the  turbulent 
shear  flow  and  a mechanically-generated  progressive,  sinus- 
oidal water  wave  are  described.  The  data  acquisition 

and  reduction  schemes  are  also  described. 

c.  in  Chapter  4,  the  experimental  results  concerning  the 
amplitude  and  phase  distributions  of  wave-induced  velo- 
cities and  turbulent  Reynolds  stresses  across  the  turbulent 
boundary  layer  are  presented.  The  importance  of  wave- 
induced  turbulent  Reynolds  stresses  on  the  momentum  trans- 
fer from  wind  to  wave  is  verified  based  on  the  experimen- 
tal data. 

d.  In  Chapter  5,  the  turbulent  models  of  Norris  and  Reynolds 
(1975)  , as  adapted  for  the  present  wave-turbulent  bound- 
ary layer  interaction  problem,  are  described.  Furthermore, 
the  comparison  between  the  results  of  model  computations 
and  the  experimental  data  is  presented. 


e.  In  Chapter  6,  we  summarize  our  work  and  present  the  con- 
clusions and  recommendations. 


T 


2.  THEORETICAL  BACKGROUND 


A periodic  disturbance  in  a turbulent  shear  flow  will  intro- 
duce additional  periodic  velocities  and  pressure  in  the  flowfield. 
In  this  study,  we  consider  turbulent  air  flow  and  a mechanically- 
generated,  single  frequency , progressive  water  wave  at  the  air 
water  interface,  as  a coupled  dynamic  system.  Our  interest  is  to 
study  the  wave-induced  motions  in  the  turbulent  shear  flow.  Two 
distinct  averaging  procedures  are  used,  as  defined  by  Hussain  and 
Reynolds  (1970),  to  sort  out  the  statistical  contribution  of  the 
periodic  disturbance  from  the  background  turbulent  fluctuations. 

The  averaging  procedures  enable  us  to  decompose  the  air  flow  into 
overall  mean,  wave-induced,  and  turbulent  components.  Applica- 
tion of  these  averaging  procedures  to  the  momentum  and  continuity 
equations  yields  the  resulting  equations  governing  the  wave-induced 
quantities.  The  appearance  of  the  wave-induced  turbulent  Reynolds 
stresses  in  the  equation  and  their  relation  to  the  momentum  trans- 
fer from  wind  to  wave  is  also  examined. 


2 . 1 Decomposition  and  Averages 

The  instantaneous  signal  f(x,t)  in  a turbulent  air  stream 
disturbed  by  a small  amplitude  progressive  water  wave  can  be 
decomposed  into  three  different  components  (see  Figure  2.1); 


f(x,t)  = f(x)  + f'(x,t)  + f ( x , t ) 


(2.1) 
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Here  f (x)  is  the  time- independent  mean  component,  f(x,t)  is  the 
wave-induced  quantity,  and  f'(x,t)  is  the  background  turbulence. 
In  order  to  single  out  tiie  wave-induced  quantities,  two  averaging 
procedures  are  defined.  The  time  average  which  yields  the  mean 
value  of  f(x)  is  defined  as 


1 ' 

f (x)  = lint  y f (x, t) 

™ Jo 


dt 


(2.2) 


The  phase  average  is  defined  as  an  average  of  the  signal  at  a 
particular  phase  of  the  wave.  Mathematically  this  can  be  ex- 
pressed as 


1 N 

f(x,t)'  = lim  — Y.  f (x,  t + n i ) 
— N 

n=o 

N-*  00 


wherex  is  the  period  of  the  wave-induced  disturbance.  The  back- 
ground turbulence  is  assumed  to  be  random  and  to  make  no  contri- 
bution to  a phase  average.  Hence,  the  phase  average  of  f(x,t) 
contains  only  the  mean  and  wave-induced  parts,  viz., 


<f(x,t)>  = f (x)  + f ( x , t ) 


(2.3) 


The  wave-induced  part  can  be  obtained  by  subtracting  the  mean 
quantity  from  the  phase-averaged  quantity;  hence. 


f(x,t)  = <f(x,t)>  - f(x) 


(2.3a) 


The  wave-induced  quantity  f(x,t)  describes  the  interaction  be- 
tween the  turbulent  shear  flow  and  the  progressive  water  wave. 
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The  equations  governing  the  wave-induced  velocities  u^  and  pres- 
sure p are  derived  in  the  following  section. 


2.2  The  Governing  Equations 

The  equations  governing  the  incompressible,  constant  property, 
two-dimensional  turbulent  boundary  layer  air  flow  in  the  wind-wave 
channel  (as  shown  in  Figure  3.1)  are  the  continuity  and  momentum 
equations.  The  dimensionless  form  of  the  equations  based  on  the 
free-stream  velocity  and  the  channel  half  height  D (air  portion 
only)  are 


9u.  9u. 

i , i 

9t  Uj  9x. 

J 


Re  9x.9x. 
J 3 


where  Re  = D/  v is  the  Reynolds  number  and  v is  the  kinematic 
viscosity.  The  velocity  and  pressure  fields  in  the  air  can  be 


decomposed  into  three  components  as 


u.  = u . + u . + u ' 
1111 


(2.6) 


p = p + p + p' 


(2.7) 


We  substitute  Equations  (2.6)  and  (2.7)  for  u.  and  p in  Equa- 
tions (2.4)  and  (2.5).  Taking  the  phase  average  and  time  average 
of  the  resulting  equations  and  subtracting  the  time-averaged  , 


E2KK&*  v & £F'%- 


-r  '-rr*  . 


equations  from  the  phase-averaged  equations  yields  the  equations 
governing  the  wave-induced  quantities  (Hussain  and  Reynolds 
1970). 


3u . 

— — — — 0 (2.8) 
dX  . 

i 


_ - 2' 

3u.  3u.  . 3u.  3p  l 3 u. 

i + — l , i = - , L_ 

3t  uj  3x.  Uj  3x.  3x.  Re  3x.3x. 

J 3 J 1 i J J 


+ (u.u.  - u.u.)  - -r— — (<uTur>  -uTuO  (2.9) 

"X.  1 J l J 3Xj  l j l j 


In  a small-amplitude,  linear  analysis  the  term 


a 

3x . 
J 


[u.u. 
1 J 


u.u.) 
1 J 


is  quadratic  in  the  induced  quantities  and  can  be  neglected. 
The  term 


<ufu'>  - uTuT  = r . . 

13  13  13 

is  the  difference  between  the  phase  averaged  and  time  averaged 
Reynolds  stresses  and  represents  the  oscillating  components  of 
the  turbulent  Reynolds  stresses  caused  by  the  passage  of  waves. 

Having  neglected  the  nonlinear  term  in  Equation  (2.9),  we 
produced  in  Equations  (2.8)  and  (2.9)  a linear  system  in  the  in- 
duced quanties  u.  and  p.  The  r..  terms  are  unknown:  hence,  a 
i 13 

closure  assumption  must  be  made.  Both  inviscid  and  quasi- 
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laminar  models  solve  the  closure  problem  by  assuming  that  the 
travelling  waves  do  not  influence  the  turbulent  structure  and, 

hence,  that  r.^  is  zero.  Miles  inviscid  model  neglects  the 

2~ 

• i 3 u . . 

viscous  term  1 , also. 

Re  3x.3x. 

J J 


To  seek  a closure  condition  for  r..,  one  could  derive  an 

ij 

equation  for  r„  by  employing  the  equations  for  the  background 
turbulence  (Hussain  and  Reynolds,  1970),  namely. 


3u(  3u^  _ 3u(  3u.  3u. 

1 . — l,  l . l , l 

-r— - + u . -r — + u . - — + u . - — + u . - — 

3t  i 3x.  1 3x.  1 3x.  1 dx. 

J J 3 J J 


2 

3p'  j 

+ r- — + — (<uTu(>  - U (u  .') 

3x.  Re  3x.3x.  3x.  l j l y 

i J 3 J 


(2.10) 


: 


When  the  governing  equation  for  r are  derived,  more  unknowns 
appear  in  the  equations.  Consequently,  the  closure  problem  is 
pushed  to  even  higher  orders.  We  are  forced,  then,  to  resort  to 


a semi-empirical  turbulent  closure  models  by  relating  r_  to  the 
other  perturbation  quantities.  The  eddy  viscosity  model  considers 
that  the  effect  of  the  turbulent  mixing  is  similar  to  the  effect 
of  viscosity  in  that  r is  proportional  to  the  wave-induced 
velocity  field  through  an  eddy  viscosity.  For  "higher  order"  tur- 
bulent energy  closures  (see  Reynolds,  1970),  one  can  employ  an 

~2  '9 

equation  for  the  wave-induced  turbulent  kinetic  energy  q = u.v 
in  conjunction  with  constitutive  equations  relating  r . ^ to 
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The  wave-induced  turbulent  kinetic  energy  equation  can  be 


derived  through  procedures  similar  to  those  used  in  obtaining  the 
r equation.  The  eddy  viscosity  and  turbulent  energy  closure 
models  of  Norris  and  Reynolds  (1975)  adapted  for  the  present  study 
are  discussed  in  Chapter  5.  Although  the  wave-induced  equations 
are  not  closed,  the  role  played  by  r_  in  the  momentum  transfer 
from  wind  to  wave  can  be  examined  (see  the  following  section). 


2 . 3 Momentum  Transfer  from  Wind  to  Wave 

The  wave-induced  turbulent  Reynolds  stress  contribution  to 
the  momentum  transfer  from  wind  to  wave  was  evaluated  both  by 
Phillips  (1966)  and  Miles  (1967).  Phillips  (1966,  page  94)  derived 
from  the  wave-induced  momentum  equation  an  expression  for  the  mean 
rate  at  which  momentum  is  transferred  from  wind  to  wave,  say  F 
per  unit  area,  i.e., 


- 3n  3r!2 

F = - pa(uv)n  - Pa(rn--  - n 3^-  )n 


(2.11) 


where  the  subscript  n indicates  these  quantities  are  evaluated  at  the 

air  water  interface  and  n is  the  displacement  from  the  mean  of  a 

periodic,  surface  wave  propagating  in  the  horizontal  x-direction. 

The  wave-induced  Reynolds  stress  - p (uv)  can  be  expressed  in 

3 n 

terms  of  wave-induced  vorticity  fl  as 


P (uv)  = p I vtt  dy 
a n a I 

> n 


(2.12) 


f.  *#  ■*•*»*■  , 


(t- 


where 


, _ 3u  3v 

3y  3x 


and  can  be  related  to  r _ through  the  equation  for  the  wave- 
induced  vorticity  (Phillips  1966,  page  97) 


2_ 

. 3ft  . d u 

(u-c)  ix  + 

dy 


3 3,  N , . 3 

v = — (r  - r ) + ( — 

3x  9y  22  il  3x2 


2 2 
3 


3y 


- ) r 

2 1 2 


(2.13) 


through  use  of  the  following  procedures: 

„ , , ...  ~ i'i  i(Kx-wt+;  1 

a.  Substitute  the  travelling  wave  expressions(f= | f | e f') 

for  ft,  v and  r„  into  Equation  (2.13)  and  express  ft  as 

a function  of  r..,  v and  u. 
ij 

b.  Substitute  the  result  of  a.  into  Equation  (2.12). 

c.  After  simplifying,  express  the  wave-induced  Reynolds  stress 
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Here  p is  the  air  density,  K is  the  wave  number  and  the  expres- 
a 

sion  for  vco  in  Equation  (2.14)  is 


— = 1 
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(2.15) 


where  0 and  9 . are  the  phase  lags  between  v and  n and  r and  n, 

2 i.1 

respectively.  Let 


d u 2 
2 V 

d_Z 


(2.16) 
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3n  STj 

and  F =-  P (r  n - ' ) 

w2  a 11  3x  ay  n 


(2.18) 


Thus,  equation  (2.11)  can  be  expressed  as 


F = F + F + F 

c wl  w2 


The  momentum  transfer  F given  by  (2.16)  is  the  contribution  eval- 

c 

uated  at  the  critical  height  (where  the  mean  velocity  is  equal 
to  the  wave  celerity).  The  Fc  transfer  term  is  identical  in  form 
with  that  for  Miles'  inviscid  model,  but  has  the  significant  dif- 
ference that  v (the  wave  induced  air  velocity  in  the  veritical  y 
direction)  depends  implicitly  on  r . . . Equations  (2.17)  and  (2.18) 
represent  additional  sources  of  mean  momentum  transfer  to  the 

wave.  F , is  the  contribution  due  to  r.,  at  all  layers  above 
wl  1J 

the  air-water  interface  except  at  the  critical  height.  F ^ is 
the  contribution  due  to  r..  at  the  air-water  interface. 

ij 

Since  r and  r are  sinusoidal  functions  of  x,  the  travelling 

11  12 


i(Kx-o)t+e  ) 

wave  expressions  r = r e 11  and 

11  'll1 


, . i(Kx-wt+e  ) , , . „ _ , » i □ T. 

r = r e \ ? can  be  used  in  Equation  2.18.  The  re- 

12  1 1 ?' 


resulting  equation  is 
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(2.19) 


For  the  invlscid  and  quasilaminar  models,  the  turbulence  structure 

is  assumed  not  to  be  affected  by  the  wave  and,  thus,  there  are 

no  wave-induced  components  of  the  turbulent  Reynolds  stresses. 

Therefore,  both  F , and  F „ are  zero.  The  importance  of  r..  in 
wl  w2  tj 

the  momentum  transfer  from  wind  to  wave  can  be  evaluated  ac- 
cording to  Equations  (2.14),  (2.15)  and  (2.19),  by  use  of  experi- 
mentally measured  values  of  the  wave-induced  velocity  and  turbulent 
Reynolds  stresses.  The  experimental  program  conducted  to  study 
the  interaction  between  a turbulent  shear  flow  and  a progres- 
sive water  wave  is  described  in  Chapter  3.  The  methods  used  to 
obtain  both  the  magnitude  and  phase  of  wave-induced  velocity  and 

r..  are  also  discussed  there, 
ij 
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3.  EXPERIMENTAL  METHODS  AND  DATA  REDUCTION 


In  this  chapter,  the  experimental  facilities  and  instrumen- 
tation used  in  this  study  are  described.  The  velocity  and  wave 
height  signals  were  measured  by  hot  film  sensors  and  a capacitance 
wave  height  gauge,  respectively.  The  hot  film  probe  was  mounted 
on  a vertical  traversing  mechanism  to  allow  measurement  throughout 
the  whole  turbulent  boundary  layer.  A wave  following  system  was 
also  used  in  this  study  to  supply  data  in  close  proximity  to  the 
air-water  interface.  Data  were  processed  and  analyzed  on  a digi- 
tal data  acquisition  and  reduction  system.  The  wave-induced 
perturbations  in  the  turbulent  air  stream  were  singled  out  by 
using  both  cross  correlation  and  cross  spectrum  techniques. 

3 . 1 Experimental  Facilities 

3.1.1  Wind  Wave  Channel 

The  experiment  was  conducted  in  the  Stanford  Wind,  Water- 
Wave  Research  Facility  as  shown  in  Figure  3.1.  A detailed  descrip- 
tion of  the  facility  was  given  in  Hsu  (1965).  The  overall  length 
of  the  channel  is  approximately  40m.  The  test  section  is  approxi- 
mately 20m  long,  lm  wide,  and  2m  high.  Water  depth  in  the  channel 
is  variable.  For  this  study  the  channel  was  filled  with  water  to 
a lm  depth  and  the  upper  half  of  the  channel  served  as  the  air  flow 
section.  A horizontal-displacement-type  mechanical  wave-generator 
capable  of  generating  selected  wave  forms  is  located  at  the  up- 
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stream  end  of  the  channel.  In  the  present  study,  a sinusoidal  sig- 
nal of  frequency  1 Hz  from  a function  generator  with  a frequency 
resolution  of  + 1%  formed  the  input  to  the  closed-loop  servo-control 
system  for  the  generation  of  a progressive  wave  with  an  amplitude 
of  2.54  cm.  A basket  consisting  of  stainless  steel  turnings  is 
placed  behind  the  wave  plate  to  minimize  wave  reflection.  A beach 
consisting  of  stainless  steel  turnings  is  placed  at  the  downstream 
end  of  the  channel  for  the  same  purpose.  The  wave  reflection  coef- 
ficient for  waves  at  a frequency  of  1 Hz  was  found  to  be  about  5 
percent  (Bole  and  Hsu,  1967).  Air  flow  is  produced  by  drawing  air 
through  the  test  section  with  a suction  fan  at  downstream  end  of 
the  channel.  A stainless  steel  honeycomb  is  placed  before  the  fan 
to  eliminate  vortex  motion  of  the  air  flow  caused  by  circular  motion 
of  the  fan.  The  air  inlet  is  positioned  approximately  4m  down- 
stream of  the  wave  plate  to  allow  full  establishment  of  the  mechani- 
cal waves  prior  to  wind  action.  The  air  moves  through  fiberglass 
filters,  a honeycomb  and  3 wire  mesh  screens  before  entering  into 
the  test  section.  A flat,  smooth  aluminium  plate  approximately  2m 
long  was  installed  at  downstream  of  the  air  inlet  to  provide  a smooth 
transition  to  the  water  surface.  The  cross-sectional  distribution 
of  air  flow  in  the  channel  reported  by  Hsu  (1965)  indicates  that 
flow  is  two  dimentional  for  substantially  the  entire  test  section. 

The  location  of  the  measurement  station  was  selected  at  approxi- 
mately 13m  from  the  end  of  the  aluminum  plate  in  front  of  the  air 
inlet. 
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3.1.2  Probe  Carriage  and  Traversing  Mechanism 

The  probe  carriage  and  traversing  mechanism  permit  re- 
mote contral  of  probe  location  in  the  longitudinal,  transverse, 
and  vertical  directions  within  the  wind-wave  facility.  The  reso- 
lution in  measuring  the  vertical  distances  is  0.0254cm.  Details 
of  the  motorized  carriage  and  traversing  mechanism  are  given  by 
Chambers,  et  al  (1970).  A sketch  of  the  probe  carriage  and  sup- 
port is  shown  in  Figure  3.2.  The  fixed  frame  measurements  were 
conducted  with  the  velocity  probes  mounted  on  the  traversing  me- 
chanism. Data  were  taken  at  ^ ixed  heights  between  3.2cm  and 
34cm  above  the  mean  water  surface. 


sembly,  and  a stainless  steel  tube  (used  to  attach  the  sensors) 


that  follows  the  wave.  The  bulk  of  the  wave-follower  was  installed 
outside  the  wind  wave  channel  and  on  top  of  an  aluminium  struc- 


i ture,  as  shown  in  Figure  3.3.  Only  the  supporting  aluminium 

> 

S channel  extended  into  the  channel.  The  electrical  controls  of 

t 

I 

the  wave  follower  system  include  (i)  a position  potentiometer  for 


which  provides  meters Avoltage  scaling  in  according  to  the  wave 


height  gauge  calibration.  The  electric  and  mechanical  devices  are 
coupled  in  the  wave  follower  system  through  a negative  feedback 
control  system.  This  system  is  capable  of  a maximum  oscillation  of 
12.7cm.  The  signal  from  the  capacitance  wave  height  gauge  provides 
the  input  to  the  mechanical  system,  changes  in  water  surface  eleva- 
tion creating  an  error  voltage  between  the  position  potentiometer 
and  the  input  of  wave  gauge.  This  error  is  amplified  and  used 
to  correct  the  system  location,  thereby  restoring  correspondence 
between  input  and  output  signals.  Details  of  the  mechanical  and 
electrical  control  systems  were  described  by  Yu,  et  al  (1971, 

1973) . A tachometer  was  attached  to  the  back  end  of  the  shaft  of 
the  motor  to  indicate  its  speed  of  rotation.  The  rotation  speed 
of  the  motor  corresponds  directly  to  the  oscillating  velocity  of 
the  wave  follower.  A calibrated  relationship  of  the  voltage 
output  from  the  tachometer  and  the  oscillating  velocity  of  the 
wave  follower  was  used  to  subtract  the  wave-follower-induced  ve- 
locity from  the  hot-f iim-sensed  vertical  velocity  component  sig- 
nal. In  this  study,  the  calibrated  result  was  1 volt  per  0.305 
m/s.  The  phase  lag  of  the  wave  follower  feedback  control  system 
at  wave  frequency  1 Hz  is  2°  (calibrated  by  cross  correlating 
between  wave  and  wave  follower  signals).  The  phase  lag  was  com- 
pensated by  placing  the  wave  gauge  wire  0.864cm  (equivalent  dis- 
tance of  2°  phase  angle  for  1 Hz  water  wave)  upstream  from  the 


other  measuring  probes.  After  compensation  the  wave  follower 


T 


motion  was  in  phase  with  the  wave  and  the  follower  maintained 
at  a selected  distance  from  the  water  surface.  The  data  obtained 
in  the  wave-following  frame  of  reference  were  taken  at  8 different 
heights  between  1.02  cm  and  5.33  cm  above  the  instantaneous  water 
surface. 


3.2  Instruments 


3.2.1  Hot  Film  Probe 

The  horizontal  and  vertical  turbulent  velocity  com- 
ponents were  measured  using  quartz-coated  cross  hot-film  sensors 
of  diameter  0.05  mm  and  length  1.016  mm.  The  hot  films  were 
operated  in  a constant  temperature  mode  and  driven  by  Thermo- 
Systems  Inc.  Model  1010  anemometers.  The  film  acts  as  one  leg 
of  the  bridge  circuit.  The  effective  cooling  velocity  (velocity 
normal  to  the  hot  film)  causes  the  film  resistance  to  change.  The 
corresponding  voltage  output  due  to  the  resistance  change  is  am- 
plified and  used  as  a feedback  signal  in  maintaining  the  bridge 
balance.  During  the  experiment,  the  films  were  operated  at  a 
sufficiently  high  over  heat  ratio  (1.6)  to  minimize  the  possibi- 
lity of  temperature  contamination  of  the  velocity  signals.  The 
experimental  correlation  applicable  to  hot  films  is 

E2  = A + Bu" 

ef  f 


where  E 

is  the 

voltage 

drop 

across 

the  wire,  U rr  is  the 
ef  f 

cool ing 

velocity 

norma  1 

to  the 

film. 

and  A, 

B and  n are  constants 

obtained 
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bv  least  square  curve  fit. 

Calibration  of  the  hot  film  sensors  was  made  in  the  core  flow 
region  of  the  wind-wave  channel,  immediately  before  and  after  each 
run.  The  mean  velocity  measured  by  the  hot-film  was  also  cross- 
checked with  that  measured  by  a pitot-static  probe,  the  difference 
between  the  two  results  being  within  3 percent  of  the  mean  velocity. 
In  order  to  accurately  calculate  the  effective  cooling  velocity, 
the  orientation  of  the  hot  films  relative  to  the  probe  centerline 
were  measured  with  an  optical  comparator  and  the  angle  between 
the  probe  axis  and  the  mean  flow  direction  was  obtained  by  suc- 
cessive calibrations  with  the  probe  rotated  180°  about  the  probe 
axis.  The  angle  between  the  hot  film  and  probe  centerline  was 
found  to  be  within  + 7°  from  the  nominal  angle  of  45°  for  the  hot 
films  used  in  this  study.  The  angle  between  the  probe  axis  and 
the  mean  flow  direction  was  found  to  be  less  than  3°  after  a care- 
ful experimental  setup.  Figure  3.4  presents  typical  cross  hot 
film  calibration  curves. 

3.2.2  Pitot-Static  Tube 

A 0.317  cm  o.d.  United  sensor  pitot-static  tube  con- 
nected to  Pace  differential-pressure  transducer  (Model  P90D, 
full  range  of  + 2.12  cm  of  water)  and  Sanborn  Series  650  recorder 
(Model  656-1100  Carrier  Amplifier)  were  used  to  calibrate  the 
hot  film  probe  and  cross-check  the  mean  velocity  obtained  from  the 
hot  film  probe.  The  Pace  differential-pressure  transducer  was 
calibrated  with  a Combust  micromanometer  with  a resolution  of 
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+ 0.006  cm  of  fluid  of  specific  gravity  0.82.  A typical  cali- 
bration curve  is  shown  in  Figure  3.5. 


3.2.3  Wave  Height  Gauge 

A Nyclad-insulated  copper  conductor  wire  was  used  as 
the  capacitance  wave  height  gauge.  Figure  3.6  is  a sketch  of  the 
wave  gauge;  a detail  report  of  the  gauge  was  given  by  Colonell 
(1966).  The  wire  is  acting  as  one  leg  of  a Sanborn  958-1100 
capacitance  bridge  amplifier  circuit.  A change  in  water  sur- 
face elevation  causes  a change  in  the  capacitance  of  the  wire. 
Before  the  calibration,  the  wire  was  presoaked  for  several  hours 
to  achieve  a stable  condition.  Static  calibration  of  the  capa- 
citance wire  was  achieved  by  attaching  the  wire  to  the  probe 
traversing  mechanism  and  systematically  varying  the  submergence 
of  the  wire  into  still  water  in  the  channel.  Accuracy  of  the  wire 
was  within  + 1 percent  of  the  mean  wave  height.  Figure  3.7  shows 
a typical  calibration  curve  of  the  wave  height  gauge. 
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3 . 3 Data  Acquisition  and  Reduction 

3.3.1  Data  Acquisition  and  Reduction  System 

The  central  unit  of  the  data  acquisition  and  reduction 
system  is  a HP  2100A  computer  with  32K  of  16  bit  core  memory.  An 
analog  to  digital  converter  system  with  a 16  channel  multiplexer 
capable  of  a maximum  sampling  rate  of  45000  samples  per  second 
was  used  to  digitize  the  input  signals  in  a simultaneous-sample-hold 
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mode.  The  load  capacity  of  the  analog-digital  system  is  + 10  volt 
with  a 5 mv.  resolution.  A CRT  terminal  is  used  for  the  purpose 
of  interrupting  or  interrogating  the  computer.  Program  input  to 
the  controller  is  accomplished  through  a high-speed  paper  tape 
reader.  An  IBM  digital  magnetic  unit  with  compatible  9-track 
format  and  800  characters  per  inch  density  was  used  for  data  re- 
cording. The  signal  to  noise  ratio  is  approximately  1000.  An 
evaluation  report  of  this  system  was  given  by  Takeuchi  and  Mogel 
(1975). 


3.3.2  Data  Acquisition  Procedures 

The  following  were  the  steps  for  the  data  taking 

procedures : 

a.  Calibrate  the  wave  height  gauge. 

b.  Calibrate  the  hot  film  probe  and  cross-check  the  mean  velo- 
city readings  obtained  from  both  the  pitot  tube  and  the  hot 
film  probe. 

c.  Check  the  wave  tank  water  level.  A pointer  located  2.54  cm 
below  the  probes  was  made  to  contact  the  still  water  sur- 
face, thus  registering  the  initial  probe  position  above  the 
mean  water  surface. 

d.  Generate  the  mechanical  wave  of  frequency  1 Hz  and  wave  am- 
plitude 2.54  cm.  Blow  the  wind  for  about  30  minutes  to 
achieve  steady  state  in  the  channel.  Set  the  probe  as  close 
to  the  wave  crest  as  possible  (avoiding  submergence  of  the 
sensors) . 
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e.  D.C.  bias  and  amplify  the  signals  to  within  the  maximum 
dynamic  range  of  the  computer  (+  10  volts).  The  turbulent 
energy  above  500  Hz  was  analyzed  to  be  small.  Therefore, 
the  signals  were  also  low  pass  filtered  at  500  Hz  with 
Rockland  1022F  filters  to  avoid  the  high  frequency  noise 
introduced  by  the  electrical  instruments.  Monitor  both  the 
hot  film  and  wave  height  signals  with  an  oscilloscope  (Wave- 
follower-induced  vertical  velocity  signal  was  also  moni- 
tored during  the  wave-following-frame  data  taking). 

f.  Digitize  the  analog  signals  every  0.001  second  for  5 minutes 
and  store  the  data  on  the  magnetic  tape  (consists  of  600 
data  blocks;  each  block  consists  of  512  data  points  from  each 
individual  signal) . 

g.  Change  the  position  of  the  probe  and  repeat  the  data  taking 
procedures  from  step  e. 
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3.3.3  Data  Reduction  Procedures 

The  following  steps  were  used  for  the  data  reduction 

procedures : 

a.  Use  the  HP  2100A  computer  and  a FORTRAN  program,  which  includes 
the  calibration  characteristics  of  sensors,  to  convert  the 
digitized  voltages  to  the  physical  quantities  of  horizontal 
velocity,  veritical  velocity  and  wave  height  (for  the  wave- 
following-frame  data  reduction,  the  induced  vertical  ve- 
locity caused  by  the  motion  of  the  probe  was  subtracted  from 
the  computed  vertical  velocity). 
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b.  Calculate  the  mean  quantities  of  data  in  Step  a. 

c.  Subtract  the  results  of  b from  those  of  a to  get  u + u', 
v + v'and  q. 

d.  Do  cross  correlation  or  cross  spectra  between  n and  u + u' 
and  q and  v + v ' . 

e.  Do  the  computation  of  (u  + u')  * (u  + u'),  (u  + u')  x (v  + v') 
and  (v  + v')  x (v  + v'). 

f.  Do  cross  correlation  or  cross  spectra  between  n and  the  results 
of  e . 

g.  Curve  fit  the  results  of  d and  f with  a sinusoidal  curve  to 
get  the  magnitude  and  phase  of  wave-induced  velocities  and 
turbulent  Reynolds  stresses. 


The  cross  correlation  and  cross  spectrum  techniques  used  for  the 
data  reduction  in  this  study  are  discussed  in  Sections  3. 3. A and 

* ■ 

3.3.5. 

3.3.4  Cross-Correlation 

Cross-correlation  is  a technique  by  which  a signal  of 
certain  known  frequency  can  be  extracted  from  the  total  signal. 

The  definition  of  cross-correlation  is 

Rnf  (t  ) = lim  | fr  n(t)  f ( t- t) dt 
T>  J o 

• I 
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Here  R , is  the  cross  correlation  function  and  t is  a known  time 

nf 

delay.  If  n is  a sinusoidal  signal  and  f is  a signa]  containing 
a sinusoidal  component  f at  the  same  frequency  as  n with  a phase 
lag  0,  R j.  can  be  expressed  as 


n f 


cos  (to T - 9) 


(3.1) 


where  to  = 2 tt  f is  the  radian  frequency. 

In  practice  the  integration  process  is  replaced  by  a sum- 
mation, namely. 


VT  } = N ifin(ti)f(ti  -T) 


where  N is  a large  number.  Typically,  N-values  of  120,000  were 
used  in  this  study.  Figure  2.1  illustrates  the  cross  correla- 


The  magnitude  and  phase  lag  of  wave-induced  component  f 
were  obtained  by  a least  square  cosine  curve  fit  of  the  cross 
correlation  results  with  (3.1).  Typical  results  of  cross  cor- 
relation fitted  by  a cosine  curve  for  the  wave-induced  veloci- 
ties u and  v are  shown  in  Figure  3.8a.  The  corresponding  results 


for  the  induced  turbulent  Reynolds  stresses  are  shown  in  Figure 


Transform  pair. 


In  order  to  cross-check  the  results  from  correla- 


tion and  to  improve  the  calculation  speed,  cross  spectrum  calcula- 
tion with  the  Fast  Fourier  Transform  technique  was  used  also  in 
this  study  to  obtain  the  magnitude  and  phase  relation  of  wave- 
induced  quantities  relative  to  the  wave.  In  general,  the  power 
spectral  function  G j.(n)  is  complex  with  real  and  imaginary  parts 
representing  in  phase  and  90°  - out  of  phase  power  respectively; 
thus , 


Gnf(n)  = Cnf(n)  " 1 Qnf(n) 


where  C , is  the  co-spectra  and  Q c is  the  quadrature  spectra, 
n f n f 

The  magnitude  of  the  cross  correlation  function  is 


n ||f 

2 


-[ 


Vn)  + V(n) 


r 


The  phase  relation  between  n anu  f is 


, , -I  rQnf(n)l 

nf(n)  = tan 


Figure  3.9  shows  a typical  amplitude  versus  f plot  of  r , r and 

11  12 

r calculated  by  this  technique  and  using  the  measured  signals 

of  turbulent  Reynolds  stresses  and  the  waves.  The  resolution 
of  the  spectra  is  0.098  Hz.  The  results  show  that  the  magni- 
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tudes  at  higher  harmonics  of  the  wave  frequency  are  considerably 


smaller  than  that  at  the  generated  wave  frequency.  The  coherence 


between  two  signals  is  defined  as 


Gnf(n) 

Coherence  “ 

Cn(n)xGf(n) 


where  G^and  G^  are  the  auco  spectral  functions  of  n and  f respec- 
tively. The  coherence  gives  information  on  how  two  signals  cor- 


relate with  each  other.  Figure  3.10  is  the  coherence  spectrum 


of  r with  the  wave.  The  results  show  strong  coherence  at  har- 

2 2 

monies  of  the  wave  frequency.  Similar  features  were  observed  in 


the  coherence  spectra  of  r and  r with  the  wave. 

11  12 


4.  EXPERIMENTAL  RESULTS  AND  MOMENTUM 
TRANSFER  FROM  WIND  TO  WAVE 


In  this  chapter,  the  experimental  results  of  the  interaction 
between  the  1 Hz  sinusoidal  progressive  water  wave  and  the  turbu- 
lent shear  flow  are  presented.  The  magnitude  and  phase  of  wave- 
induced  velocity  and  turbulent  Reynolds  stresses  are  calculated 
both  in  the  fixed  and  wave-following  frames  of  reference. 

The  contribution  of  r..  to  the  momentum  transfer  from  wind  to 
wave  as  described  in  section  2.3  is  distributed  through  the  turbu- 
lent boundary  layer  above  the  air-water  interface.  To  verify  the 
importance  of  r ^ . in  the  momentum  transfer  from  wind  to  wave,  we 
need  the  wave-induced  velocities  and  turbulent  Reynolds  stresses 
across  the  whole  turbulent  boundary  layer.  Recall  that  the  fixed 
probe  measurements  were  conducted  to  obtain  the  data  at  20  fixed 
heights  between  3.2  cm  and  34  cm  above  the  mean  water  surface. 

The  data  obtained  in  the  wave-following  mode  was  extended  to  the 
region  as  close  to  the  wave  surface  as  possible  without  damaging 
the  sensor,  covering  the  region  from  1.02  cm  and  5.33  cm  above 
the  wave  surface.  The  comparison  between  the  results  of  calcu- 
lation of  F by  using  the  fixed  frame  data  alone  and  then  using 
w I 

both  the  fixed  and  wave-following  frame  data  shows  the  importance 

of  r..  close  to  the  air  water  interface.  F „ was  evaluated  with 
lj  w2 

the  wave-following  frame  data. 
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4 . 1 Mean  Velocity  and  Turbulence  Intensities 

Mean  velocities  over  the  progressive  water  wave  as  a function 
of  height  above  the  mean  water  surface  are  shown  in  Figure  4.1  for 
5 different  free  stream  velocities  as  tabulated  in  Table  4.1  These 
mean  velocity  profiles  were  used  for  the  model  computation  as 
described  in  Chapter  5.  The  comparison  between  the  results  measured 
with  pitot-static  tube  and  hot-film  probe  is  also  shown  in  the  same 
figure  (for  the  case  U = 1.87  m/sec).  The  difference  is  within 
3%  of  the  free  stream  velocity.  The  comparison  assures  the  ac- 
curacy and  reliability  of  the  hot  film  data-taking  and  reduction 
procedures . 


Table  4.1 

Free  Stream  Velocity  Conditions 


Critical  height  above  the  mean  water  surface 


0.883 


(Non-cr it ical  layer  case) 


1.115  12.34  cm  (Thick  critical  layer  case) 


1.34  8.57  cm  (Thick  critical  layer  case) 


1.836  cannot  be  reached  by  probe 

(Thin  critical  layer  case) 

2.283  cannot  be  reached  by  probe 

(Thin  critical  layer  case) 


As  shown  in  Figure  4.1  the  mean  velocity  profiles  deviate 
from  the  logarithmic  law  in  the  lower  portion  of  the  turbulent 
boundary  layer  compared  with  the  conventional  profile  over  a 


•r  * ".n  v 


x -im 
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rigid  flat  surface.  it  seems  reasonable  to  conclude  that  the 
departure  from  the  logarithmic  law  may  be  attributed  to  the  dif- 
ferences in  the  boundary  conditions  between  a deformable  air-water 
interface  with  wave  motions  and  a rigid  immobile  surface.  It 
has  been  discussed  previously  that  the  numerical  calculations 
based  on  various  adhoc,  closure  assumptions  are  extremely  sensi- 
tive to  the  velocity  distribution  and  the  turbulence  characteristics 
near  the  interface.  This  indicates  that  detailed  experimental  data 
of  the  flowfield  in  close  proximity  of  the  air-water  interface 
should  be  valuable  for  further  Improvement  of  wind-wave  genera- 
tion theory. 

Typical  normalized  distributions  of  turbulence  intensities 
and  Reynolds  stress  as  a function  of  non-dimensional  heights  mea- 
sured In  wave-following  and  fixed  frames  of  reference  are  shown 
in  Figures  (4.2a)  and  (4.2b),  respectively.  The  distributions  of 

| the  turbulence  quantities  in  Figure  4.2a  are  presented  as  a 

■ 

function  of  non-dimensional  heights  Kh^  over  the  progressive  water 
wave,  where  K Is  the  wave  number  and  h^  is  the  height  above  the 

instantaneous  wavy  surface.  The  magnitudes  of  the  relative  tur- 

2 2 2 2 

, bulence  intensities  u ’/U  and  v /U  ' and  of  the  turbulent  Reynolds 

uo  oo 

* — 2 

stress  - uv/11  are  consistent  with  those  obtained  in  a turbulent 

CO 

\ boundary  layer  over  a rigid  flat  surface  (Klebanoff,  1954).  The 

j 

* distributions  of  the  turbulence  quantities  measured  in  a fixed 

» 

frame  of  reference  as  shown  in  Figure  4.2h  are  plotted 

A 

1 


as  a function  of  non-dimensional  height  above  the  mean  water 
surface.  The  relative  magnitude  of  the  u-component  in  the  fixed 
coordinates  is  larger  than  that  in  the  wave-following  coordinate 
whereas  the  relative  magnitudes  of  the  and  - uv  components  are 

smaller.  Plausible  reasons  responsible  for  the  differences  in 
magnitude  measured  in  different  frame  of  references  are  offered 
as  follows:  First,  if  the  mean  velocity  profile  is  assumed  to 

follow  the  undulating  wavy  surface  for  a small  amplitude  wave, 
the  fixed  probe  would  sample  different  portions  of  the  local 


the  correlation  and  cross  spectrum  techniques  are  in  good  agree- 
ment. Since  the  Fast  Fourier  Transform  subroutine  of  the  cross 
spectrum  technique  is  four  times  faster  than  the  correlation  tech- 
nique, the  cross  spectrum  technique  was  adopted  for  most  of  the 
computation  in  this  study. 

When  U /C  is  equal  to  0.883,  there  is  no  critical  height 

IX) 

in  the  flowfield.  As  shown  in  Figure  4.4,  the  phase  lags  of  u 
and  \ are  close  to  180°  and  270°  respectively;  the  phase  varia- 
tions of  u and  v across  the  turbulent  boundary  layer  are  small. 

U /C  = 1.115  corresponds  to  the  case  when  the  critical  height 
is  at  12.38  cm  above  the  mean  water  surface.  As  shown  in  Figure 
4.6,  there  is  a 170°  phase  shift  of  u near  the  critical  height. 

As  a guide  for  the  experimental  phase  lag  data  of  u and  v in  the 
region  close  to  the  air  water  interface,  the  theoretical  phase 
lags  of  u and  v at  the  mean  water  surface  were  obtained  from  the 
kinematic  boundary  conditions  for  a small  amplitude  wave  as  fol- 
lows. The  kinematic  boundary  conditions  at  the  air-water  inter- 
face are 


u | = u + CKn  + 0(nO  = u [ + n n + o(n  ) (4.1) 

'y=n  d y=0  ay'y=0 


v | = + u + 0(rr)  = v|  n + | n n + 0(n?)  (4.2) 

y=n  at  ax  'y=0  ay  'y=0 


Here  is  the  drift  velocity  at  the  air-water  interface,  and  the 
velocities  at  the  air-water  interface  are  expressed  in  terms  of 
the  velocities  at  the  mean  water  surface  by  using  the  Taylor 
series  expansion.  Using  the  decomposition  of  F.quation  (2.1), 
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phase  averaging  Equations  (4.1)  and  (4.2),  time  averaging  Equa- 
tions (4.!)  and  (4.2)  and  subtracting  the  time  averaged  quanti- 
ties from  the  phase  averaged  quantities  result  in 


(4.3) 


_ _3n  — i _3n 

y=0  3t  11 ' y=0  3x 


(4.4) 


With  the  travelling  wave  expression  n = | n I cos (Kx-ujt)  substi- 
tuted into  Equations  (4.3)  and  (4.4)  we  have 


In  this  study  a)  = 2 tt  (frequency  = 1 Hz)  and  K = 2 tt  /C ; thus.  Equa- 


tions (4.5)  and  (4.6)  can  be  expressed  as 


The  measured  phase  lags  of  u and  v at  the  lowest  height  of  measure- 
ment (3.18  cm  above  the  mean  water  surface)  were  compared  with  the  re- 
sults obtained  from  Equations  (A.  7)  and  (A. 8)  as  shown  in  Table  A. 2.  In 
this  experiment,  > 2 it  at  the  lowest  measurement  point  for  all  the 
experimental  runs  (based  on  the  measured  mean  velocities  at  the  lowest 
two  measurement  points  above  the  mean  water  surface).  Experimental  data 
also  show  that  u < C for  the  cases  of  U /C  = 0.883,  U /C  = 1.115  and 

CO  oo 

U / C = 3. A,  while  u > C for  the  cases  of  U /C  = 2.28,  at  the  lowest 

OO  oo 

measurement  point. 

Table  A. 2 

Comparison  of  Estimated  and  Measured  (Fixed  Probe  Data)  Phase 
Lags  of  u and  v at  the  Lowest  Measurement  Point 


u 

<30 

(m/sec) 

U /C 

OO 

u„(0  ) 
e r 

.Equation. 
( A.  7 ' 

vv 

(Measured) 

vv 

.Equation. 
( A. 8 ’ 

ve 

(Measured) 

1.381 

0.883 

180° 

18A . 2A° 

270° 

273.77° 

1.  7AA 

1.115 

180° 

187.2° 

2 70° 

277.99° 

2.09b 

1.3A 

180° 

166.76° 

270° 

273. A6° 

2.872 

1.8A 

180° 

130. 79° 

90° 

2A.A5° 

3.  571 

2.28 

180° 

128.68° 

90° 

5A.A1° 

As 

shown  in 

Table  A. 2, 

the  discrepancy  between  the  experimental 

results 

and  the 

results  obtained  from  Equations  (A. 7) 

and  (A. 8)  in- 

creased 

as  the  wind  velocity 

increased . 

A plausible 

reason  for  the 

consistent  deviations  could  be  that  when  the  wind  velocity  is  in- 
creased the  shear  velocity  is  increased,  therefore  the  lowest  measure- 
ment point  in  y+  coordinates  is  raised  until  Equations  (A. 7)  and 
(A. 8)  are  no  longer  valid. 


As  shown  in  Figures  4.3,  4.5,  4.7,  4.9  and  4.11  the  magnitude 
of  u is  much  larger  than  that  of  v close  to  the  air-water  inter- 
face and  the  magnitude  difference  between  u and  v decreases 
as  distance  increases  from  the  interface.  This  may  be  due  to 
the  fact  that  the  fixed  probe  traversed  different  portions  of  the 
local  mean  velocity  profile  and  the  magnitude  of  u is  also  contri- 
buted by  the  velocity  gradient  of  the  portion  of  the  local  mean 
velocity  profile  traversed  by  the  sensor;  therefore,  as  the  mean 
velocity  gradient  is  decreasing  with  Kh  so  does  the  magnitude 
difference  between  u and  v.  As  shown  in  Figure  4.7,  the  measured 
magnitude  of  u is  decreasing  with  Kh , while  the  magnitude  of  v is 
increasing  with  Kh  near  the  interface.  The  reason  could  be  due 
to  the  fact  that  the  magnitude  of  v is  a function  of  both  the  de- 
cay factor  e ^ and  the  difference  between  the  local  mean  velo- 
city and  wave  celerity  (u  - C)  (See  Phillips  1966). 

The  measured  values  of  the  amplitude  and  phase  of  r , r 

11  12 

and  r as  a function  of  non-dimensional  height  above  the  mean 

22 

water  surface  are  shown  in  Figures  4.13  through  4.22.  As  shown 
in  Figures  4.13  (a,  b,  c),  4.15  (a,  b,  c) , 4.17,  4.19  and  4.21 

(a,  b,  c) , the  magnitude  of  j r j is  an  order  of  magnitude  larger 

1 1 

than  those  of  | r | and  ]r  | and  the  magnitude  of  |r  | is  the 

1 2 22  22 

smallest  close  to  the  air-water  interface.  The  distributions  of 

|r  | across  the  turbulent  boundary  layer  have  the  same  decreasing 

11 

trend  as  a function  of  Kh  for  all  cases,  i.e.,  non-critical  layer. 

Figure  4.13a;  thin  critical  layer.  Figures  4.19  and  4.21a;  and 

thick  critical  layer.  Figures  4.15a  and  4.17. 
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The  magnitude  and  phase  angle  relations  between  the  different 
r components  at  the  air-water  interface  may  be  estimated  as  fol- 

ij 

lows.  By  assuming  that  the  streamline  follows  the  wavy  surface, 
i . e . , that 


3n  _ v 
3x  (u-C) 


(v  + v + v ') 

(u  - C + u + u ') 


3n 


or  v + v + v'  = (u-C+u+  u')  (4.9) 

we  have  a basic  boundary  condition.  Now, 
a-  Multiply  Equation  (4.9)  by  v' . 

b.  Phase  average  the  results  of  a. 

c.  Time  average  the  results  of  a. 

d.  Subtract  the  result  of  c from  that  of  b to  obtain 


2 n | n I 


22 


1 2 


and  0 =90+0 

22  1 2 


(4.10) 


In  this  study  |n|  = 2.54  cm  and  X = 1.56  m;  therefore. 


r = 0.102  r 

22  1 2 


(4.11) 


I < 

1 | 

1 » 

I 

) 


e.  Multiply  Equation  (4.9)  by  u'. 

f.  Phase  average  the  result  of  e. 

g.  Time  average  the  result  of  e. 

h.  Subtract  the  result  of  g from  that  of  f to  obtain 


0 =90+0 

12  11 


(4.12) 
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(4.13) 


and  r | = 0.102  r 

12  11 


The  experimental  results  for  |r  | at  the  lowest  measure- 
ment point  as  tabulated  in  Appendix  II  indicate  comparable  order 
of  magnitude  among  the  three  components  of  |r„|  for  experimental 
and  theoretical  estimates,  particularly  for  the  high  wind-speed 
cases,  U /C  = 1.34,  1.84  and  2.23. 

CO 

However,  the  experimentally  measured  phase  lags  of  the  r_  at 
the  lowest  measurement  point  disagree  with  the  theoretical  esti- 
mates given  by  Equations  (4.10)  and  (4.12).  The  inconsistancy  could 
be  due  to  a sharp  gradient  of  the  phase  angle  close  to  the  water 
surface,  so  that  the  measured  data  at  the  lowest  measurement  point 
does  not  represent  the  actual  phase  angle  at  the  air-water  inter- 
face . 


4.2.2  Wave-Fol lowing  Measurements 

Figures  4.23  and  4.24  show  the  magnitudes  and  phase- 
lags  of  wave-induced  velocity  components  u and  v as  a function  of 
non-dimensional  height  Kh^.  above  the  instantaneous  wave  surface 
for  five  different  mean  free-stream  velocities.  To  establish  the  valid- 
ity of  the  experimental  data,  the  experimental  results  of  v at  the 
lowest  measurement  point  (hf=1.143  cm)  were  compared  with  the  esti- 
mate of  v obtained  through  the  kinematic  boundary  condition  on  the 
wave  surface. 

Considering  a sinusoidal  wave  n - ( n|  cos  uit , we  find  that  the 
kinematic  boundary  condition  for  v at  the  wave  surface  is 
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= — = - ai|n|sinu>t  = tii|  ri|  cos  (oat  - 270°) 

The  theoretical  phase  lag  of  v is  270°  and  the  magnitude  of  v is 
u)|n|  at  the  wave  surface.  In  this  study,  the  wave  amplitude  is 
2.54  cm  and  the  wave  frequency  f is  1 Hz.  Therefore,  the  theore- 
tical value  of  | v j is  16.0  cm/sec  at  the  wave  surface.  The  mea- 
sured magnitude  and  phase  of  v at  the  lowest  measurement  point  is 
tabulated  in  Table  4.3  for  comparison  with  the  theoretical  esti- 
mate of  v at  the  wave  surface. 

Table  4.3 

Measured  Values  (Wave-Following  Data)  of  Wave-Induced  Vertical 
Velocity  at  the  Lowest  Measurement  Point 


u /c 

CO 

(y+)  . 

min 

M 

V6 

1.11 

55.9 

13.045 

268.97° 

1.  34 

72.9 

10. 394 

271 .46° 

1.49 

89.1 

8.352 

276. 38° 

1.84 

128.0 

5.09 

291.77° 

2.28 

164.0 

2.652 

303.99° 
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Although  the  distance  above  the  instantaneous  wavy  surface  is 
cLoselv  the  same  for  all  free  stream  velocities,  the  relative  posi- 
tion of  the  probe  is  farther  away  from  the  surface  in  y+-  coordi- 
nates at  high  wind  speeds  because  of  the  increase  in  u* . In  other 
words,  the  probe  was  not  able  to  measure  the  "true"  wave-induced 
quantities  at  high  wind  speeds.  It  can  be  seen  from  Table  4.3  that 
the  discrepancy  between  the  experimental  results  and  the  theore- 
tical estimates  increased  as  the  wind  velocity  increased. 

As  shown  in  Figure  4.23,  the  magnitude  of  u decreases  with 
Khj.  for  all  wind  speeds.  However,  |v|  decreases  with  Kh^  for  the 
low  wind  cases  (U  /C  = 1.11,  1.34,  and  1.49)  and  increases  for  the 

oo 

high  wind  cases  (l!  /G  = 1.84  and  2.28).  The  differences  in  | v | 
variation  with  respect  to  Kh  observed  for  the  high  and  low  wind 
speeds  may  be  due  to  the  fact  that  | v | i s dependent  on  both  the 
disturbance  decay  rate  in  the  vertical  direction  and  the  local 
value  of  (u  - C) . As  shown  in  Figure  4.24,  the  measured  phase 
angles  of  v are  close  to  the  theoretical  boundary  condition  value 
of  270°  phase  lag. 

In  order  to  show  the  differences  in  magnitude  of  the  measured 
u and  v components  in  both  the  fixed  and  wave-following  frames 
of  reference,  the  | u|  and  | v J at  the  lowest  measuring  position  as 
well  as  at  a fixed  non-dimensional  height,  Kh  = Kh^  = 0.22  were 
tabulated  in  Table  4.4. 
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Tab  1 e 4.4 


Comparison  of  u and  v Data  Both  in  the  Fixed  and  Wave  Following  Frames 


l.il  1.34  1.84  2.28 


F 

(Lowest  Point) 

7.03 

0.58 

6.45 

0.23 

3.53 

0.82 

2.77 

1 . 1 

W-F 

(Lowest  Point) 

5.04 

7.43 

3.95 

5.11 

8.31 

1.68 

7.0 

0.  54 

F 

(Kh  = 0.22) 

3.96 

0.44 

4.78 

0.23 

0.86 

0.82 

1.23 

1 . 1 

W-F 

(Khf  = 0.22) 

2.17 

6.65 

1.61 

4.10 

4.68 

2.68 

5.67 

0.79 

As  discussed  in  Section  4 . 2.  1 , assuming  that  the  mean  velocity  pro- 
file oscillates  vertically  with  the  surface  wave  motion  means 
that  the  fixed  probe  measures  different  portions  of  the  locan  mean 
velotity  profile  associated  with  wave-following  coordinates.  Con- 
sequently, the  difference  in  the  local  mean  velocity  which  is  a 
function  of  the  wave  amplitude  and  the  mean  velocity  gradient,  is 
included  in  the  measured  |u  | - component.  Since  the  velocity  gra- 
dient in  the  region  of  the  fixed-probe  decreases  as  the  wind  speed 
increases,  there  is  less  magnitude  difference  between  moving  and 
fixed  frame  values  at  high  wind  speeds  than  at  the  low  wind  speeds. 
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On  the  other  hand,  there  is  less  magnitude  difference  between  |u 


and|v|in  the  wave-following  measurements  except  at  high  wind 


speed,  IWC  = 2.23.  The  large  difference  between  |u  | and  |v  j 
at  high  wind  speed  may  suggest  "flow  separation"  where  the 


mean  velocity  profile  ceased  to  follow  the  surface  wave  motion 


as  assumed. 


Figures  4.25  and  4.27  show  the  magnitude  of  r..  as  a 


function  of  Kh^.  above  the  instantaneous  water  surface  for  five 
different  free  stream  velocities.  For  the  high  wind  cases 

(—  = 1.84  and  2.28)  |r  | is  larger  than  jr  I and  |r  I.  For 

0 11  12  22 

the  low  wind  cases,  the  orders  of  magnitude  among  r , j r 

11  12 

and  r ; are  the  same.  In  Figure  4.25  the  decay  rate  of  1 r 

22  11 

as  a function  of  Kh^  is  more  rapid  for  the  high  wind  cases 
than  it  is  for  the  low  wind  cases,  while  the  opposite  is  observed 
for  jr  j in  Figure  4.26.  In  Figure  4.27,  similar  trends  for 


r j as  a function  of  Khc  are  observed  for  different  wind  ve- 

12  f 

locities.  Figure  4.28  shows  the  phase  lags  for  r , r and 

11  12 

r . The  phase  lags  or  r..  were  found  not  sensitive  to  dif- 

22 


*i  ferent  free-stream  velocities,  and  the  variation  of  the  phase 

t 


4 . 3 Momentum  Transfer  from  Wind  _t_o  Wave 

Kendall  (1970)  studied  the  structure  of  turbulent  flow  over  a 
progressive  wavy-wall.  Because  he  did  not  make  detailed  measure- 
ments of  the  distribution  r„  (bath  in  magnitude  and  phase)  through- 
out the  turbulent  boundary  layer,  the  turbulence  contributions  of 
the  momentum  transfer  from  wind  to  waves  could  not  be  evaluated. 
However,  based  on  his  qualitative  estimation  he  stated  that  the 
induced-turbulent  Reynolds  stresses  are  probably  important  in  the 
momentum  transfer  from  wind  to  waves. 

From  the  measured  mean  velocity  profile,  wave-induced  velo- 
cities, and  turbulent  Reynolds  stresses  described  in  the  previous 
sections,  the  turbulent  contribution  to  the  momentum  transfer  from 
the  wind  to  wave  can  be  evaluated  through  the  expressions  in  Equa- 
tions (2.16)  to  (2.18)  in  Section  2.3  For  example,  from  the  experi- 
mental data  taken  at  the  lowest  measuring  point  to  the  interface, 

the  wave-induced  Reynolds  stress  at  the  interface,  - p (uv)  can 

a n 

be  estimated  and  may  be  expressed  as 


- p (uv)  = — U - — cos((jjt  + 9) 

a / 


where  6 is  the  phase  relation  between  u and  v.  The  results  of  F , 

F .,  F _ and  - p (uv)  are  tabulated  in  Table  4.5  for  5 different 
wl  w2  a n 

free-stream  conditions. 


Table  4.5 


Experimental  Results  of  Momentum  Transfer  from  Wind  to  Waves 


: 


i 


i 

* 

5 

I 

f 

l 

i 

« 

i 


11  /c 

F 

c 

(N/m' ) 

yc 

wi 

( N / m ~ ) 

w2 

(N/m2) 

-P  (uv) 

cl 

n 

1 . 1 1 

12.34 

0.13  x 10~6 

0.35  x 

IQ"6 

0.21  x 

10  4 

0.20  x 

io-4 

1.34 

8.57 

0.59  x icf5 

0.48  x 

ID’5 

0.19  x 

io~4 

0.18  x 

IQ'4 

1.49 

* 

k 

0.43  x 

10“6 

0.10  x 

IQ"4 

0.12  x 

IO"4 

1.84 

■k 

k 

0.11  x 

10-6 

0.50  x 

10“  5 

0.51  x 

IO"5 

2.28 

* 

* 

0.68  x 

10~7 

0.49  x 

IQ'5 

0.98  x 

IO"5 

* Cannot  be  measured  or  evaluated. 


For  large  values  of  Uui/C  (1.49,  1.84  and  2.28),  the  critical  height 
is  so  close  to  the  air-water  interface  that  the  critical  layer  be- 
comes inaccessible  to  the  velocity  probe.  Thus,  due  to  the  lack  of 


experimental  data  for 


d^  u 


dy 


2 > 


du 

dv 


and  v^at  the  critical  height,  F 


cannot  be  evaluated.  For  the  cases  of  U /C  = 1.11  and  1.34,  the 


d'  u 


du 


values  of  — — 7 and  — — at  the  critical  heights  were  evaluated  from 
dy  dy 


the  best  fit  of  a logarithmic  profile  to  the  measured  mean  velocity 


data.  F , is  the  contribution  due  to  r..  at  all  la' ore  above  the 

wl  1J 
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critical  height.  The  integral  of  Equation  (2.17)  was  evaluated 
from  1.02  cm  to  3.2  cm  from  the  wave  follower  measured  data  and 
from  3.2  cm  to  34  cm  from  the  fixed  probe  measured  data.  (Due 
to  lack  of  fixed  frame  data  at  U /C  = 1.49,  in  this  case  F , 

CO  W1 

was  evaluated  with  the  wave-follower  measured  data  alone.)  In 
order  to  evaluate  the  importance  of  the  wave-induced  turbulent 
Reynolds  stresses  in  the  region  close  to  the  air-water  inter- 
face, F was  also  evaluated  by  using  the  fixed-frame  data  only 
w 1 

and  the  results  are  shown  in  Table  4.6. 


Table  4.6 

Experimental  Results  of  F . (Fixed  Probe  Data) 

wi 


u Jc 

CO 

1. 

.11  1.34  1.84 

2.28 

Fwl 

0.21 

x 10"7  0.22  x io-7  0.18  x io“7 

0.19  x io-7 

The  results  of  F , in  Table  4.6  are  small  compared  with  those 
wl 

tabulated  in  Table  4.5.  One  may  conclude,  then,  the  importance  of 

r..  close  to  the  air  water  interface, 
il 

As  indicated  in  Equation  (2.14),  the  sum  of  F and  F , is 

c wl 

equal  to  the  wave-induced  Reynolds  stress  at  the  air-water  inter- 
face. As  shown  in  Table  4.5  when  U /C  is  equal  to  1.11  and  1.34 

or.)  4 

(cases  when  F can  be  evaluated  experimentally),  the  sum  of  F 
c c 

and  F , is  smaller  than  the  wave-induced  Reynolds  stresses, 
wl 
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- p^(uv)  near  the  air-water  interface.  The  difference  may  be  at- 
tributed to  the  lack  of  experimental  data  in  the  immediate  vici- 
nity of  the  air-water  interface  where  the  integral  of  F , needs 
to  be  evaluated. 

Since  the  lowest  position  of  measurement  above  the  instan- 
taneous wavy  surface  in  the  wave-following  system  is  the  same  for 
all  free  stream  velocities,  the  probe  is  relatively  farther  away 
from  the  surface  in  y coordinates  due  to  the  increase  in  u*  with 
wind  speed.  In  other  words,  the  probe  was  not  able  to  measure 
wave-induced  quantities  close  to  the  air-water  interface  at  high 

wind  speed.  Thus,  whether  the  decrease  in  F , and  - p (uv)  as 
r w2  a n 

the  wind  velocit;  increases  as  shown  in  Table  4.5  is  due  to  less 

momentum  transfer  from  wind  to  wave  is  uncertain.  Note  that  F 

w2 

has  the  same  order  of  magnitude  as  - p (uv)  , implying  that  a 

a n 

significant  contribution  to  the  total  momentum  transfer  is  de- 
rived from  the  wave-turbulence  interaction.  Clearly  r..  cannot 

tj 

be  neglected  in  the  evaluation  of  the  momentum  transfer  from  wind 
to  waves,  their  contribution  as  shown  in  Table  4.5  (U  /C  = 1.11 
and  1.34)  comprising  about  90  percent  of  the  total  transfer,  while 
the  critical  layer  contribution  is  about  10  percent. 


% tv  ’ ■>*>  <-  a 


5.  THEORETICAL  MODELING 


Both  Saeger  and  Reynolds  (1971)  and  Norris  and  Reynolds  (1975) 
studied  the  wave-induced  pressure  in  a turbulent  channel  flow  with 
a mechanically-articulated  waving  wall.  In  their  theoretical  model- 
ing, they  used  different  techniques  to  satisfy  the  time  dependent 
boundary  conditions  at  the  wavy  wall.  Saeger  and  Reynolds  solved 
the  problem  by  transformation  of  the  boundary  conditions  to  the 
undisturbed  and  mean  position  of  the  moving  wall  using  a Taylor- 
series  expansion.  Norris  and  Reynolds  solved  the  problem  by  trans- 
formation of  the  entire  problem  to  a time  dependent  coordinate 
system.  In  the  transformed  coordinate  system,  both  walls  are 
stationary  and  in  particular  the  boundary  condition  can  be  ap- 
plied exactly  at  the  moving  wall. 

Norris  and  Reynolds  (1975)  calculated  the  wave-induced  quanti- 
ties by  using  quasilaminar , eddy  viscosity,  and  turbulent  energy 
closure  models.  By  comparing  their  experimental  and  calculated  re- 
sults of  wave-induced  pressure  on  the  stationary  wall  opposite 
the  oscillating  wavy  wall,  they  concluded: 

(1)  The  coordinate  transformation  approach  is  better  than 
the  boundary  condition  transfer  approach.  Since  for  turbulent 
shear  flow  the  velocity  curvature  is  large  near  the  boundary,  ne- 
glecting higher-order  terms  in  the  Taylor-series  expan  ion  intro- 
duces errors. 

(2)  The  computational  results  in  the  transformed  coordinate 
system  checked  well  with  the  measured  stationary  wall  pressure  data 


for  upstream  running  waves  and  for  downstream  running  waves  in 
the  range  0 1 C/U^  I:  0.4.  The  eddy  viscosity  model  results  are 
slightly  better  than  the  turbulent  kinetic  energy  model  results 
compared  with  the  measured  pressure  data. 

(3)  Although  the  calculated  wave-induced  pressure  is  quite 
insensitive  to  the  wave-induced  turbulent  Reynolds  stresses 
modeling  used,  differences  in  the  calculated  wave-induced  velo- 
cities are  noticed  with  different  turbulent  models.  Thus,  the 
velocity  calculations  are  more  sensitive  to  the  assumed  models 
than  the  pressure  computations. 

Both  Benjamin  (1959)  and  Norris  and  Reynolds  (1975)  solved 
the  wave-induced  field  in  a curvilinear  coordinate  system  so  that 
the  boundary  condition  can  apply  at  the  wavy  surface.  Benjamin 
examined  the  effect  of  neglecting  viscosity  in  the  inviscid  model. 
The  results  of  his  study  indicated  that  the  viscous  effect  is  ne- 
gligible while  for  the  closure  model  r„  may  be  important.  This 
result  is  inconsistant  with  the  Norris  and  Reynolds'  conclusions. 

Norris  and  Reynolds  compared  their  stationary  wall  pressure 
data  with  a simple  inviscid  analysis.  Although  the  magnitude  of 
the  pressure  appears  to  be  inertially  dominated  for  all  wave  speeds 
(upstream  and  downstream  running  waves),  the  phase  shifts  appear 
to  deviate  from  such  a simple  inviscid  ana1 /sis  for  downstream 
running  waves,  suggesting  that  important  effects  arise  from  viscous 
and  turbulent  influences.  Inertia  dominance  appeared  also  in  their 
laminar  flow  results.  In  other  words,  the  stationary  wall  pres- 
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sure  data  cannot  be  used  as  a conclusive  test  of  the  need  for  a 

quality  of  turbulence  modeling.  Unfortunately,  no  experimental 

data  on  the  velocity  field  were  taken  in  their  study.  Therefore, 

the  present  experimental  data  on  induced  velocities  and  r..  dir- 

ij 

ectly  over  the  progressive  water  wave  can  serve  as  a more  sensitive 
test  of  their  closure  models.  The  objective  of  the  modeling  calcu- 
lation in  this  study  is  to  determine  if  either  the  coordinate 

transformation  or  the  r..  terms  are  important  in  the  wave-turbulence 

ij 

interaction  problem. 

In  this  chapter,  the  coordinate  transform,  the  eddy  viscosity 
model,  and  the  turbulent  energy  model  of  Norris  and  Reynolds  are 
described.  In  addition,  the  mean  velocity  profile,  the  eddy  vis- 
cosity and  the  mean  turbulent  kinetic  energy  across  the  turbulent 
boundary  layer  used  in  the  Norris  and  Reynolds  modeling  as  adapted 
for  the  present  wave-turbulent  boundary  layer  interaction  problem 
are  presented. 


5 . 1 Transformation  of  Coordinates 

In  the  application  of  Norris  and  Reynolds,  they  solved  the 
wave-perturbed  equations  in  the  following  new  coordinate  system: 


t*  = t 


(5.1) 


y*  = y - n(x,t)f(y) 


(5.2) 


. , , . sinh[K(2  - y)) 

W,ere  f y)  sinh(K(2  - n(x,t)  )) 

2 TT 

and  K = — is  the  wave  number. 


(5.4) 


^ . ■ t -WFV.  ^ . 4 , 


Lines  of  constant  y*  at  a given  instant  are  shown  in  Figure  5.1; 
they  are  lines  of  constant  stream  function  for  an  inviscid  flow  be- 
tween a stationary  wavy  surface  and  a flat  wall.  In  the  transformed 
coordinate  system,  the  flat  and  the  waving  surfaces  do  not  move 
relative  to  the  transformed  coordinate  system,  i.e., 

y = n(x,  t)  - y*=0 

y = 2 -*  y*=2 

Since  the  y*=0  line  is  identical  with  the  surface  of  the  propa 
gating  water  waves,  the  boundary  conditions  can  be  applied  directly 
at  the  wave  surface  in  the  x*,  y * coordinate  system.  Transferring 
the  boundary  condition  to  mean  water  surface  employing  Taylor  series 
expansion  is  not  necessary.  The  boundary  conditions  for  the  pre- 
sent wave  and  the  turbulent  boundary  layer  interaction  problem  are 
described  in  section  5.4. 


5 . 2 Eddy  Viscosi ty  Model 

The  linearized  perturbation  Equations  of  (2.8)  and  (2.9) 
applied  for  two-dimensional  boundary  layer  flow,  with  the  vertical 
mean  velocity  v and  the  x-derivative  of  u considered  to  be  small, 
are 


Cont:nuity  Equation: 


Ul  'V 

dx  <iy 


0 


(5.5) 
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x-Momentum  Equation: 
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(5.7) 


There  are  six  unknowns,  namely  u,  v,  p,  r , r and  r , in 

11  22  12 

the  above  Equations,  (5.5),  (5.6)  and  (5.7).  The  presence  of  the 
induced-turbulent  Reynolds  stresses  introduces  the  closure  problem 
for  the  set  of  equations.  The  quasilaminar  model  assumes  that 
the  passage  of  waves  produces  no  effect  on  the  background  turbu- 
lent Reynolds  stresses;  consequently,  the  r^  terms  in  the  above 
equations  are  zero.  On  the  other  hand,  if  the  interactions  between 
the  wave  and  the  turbulence  field  cannot  be  neglected,  closure 

assumptions  on  the  r..  terms  must  be  made  in  order  to  obtain  a solu- 

ij 

tion  of  the  perturbation  quantities.  The  eddy  viscosity  closure 

model  assumes  that  the  wave-induced  turbulent  Reynolds  stresses 

r are  generated  by  the  passage  of  a wave  via  the  action  of  an 

eddy  viscosity  In  addition,  the  r„  are  assumed  to  be  related 

to  the  wave-induced  strain  rate  S..  and  expressed  as 

tj 


r . . - - 2 v S . . 
ij  T ij 


(5.8) 
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where 


. 3u . 3u . 
I / i . 1 


S.  . = ~ (—  + — J) 
n 2 . 


lx . 3x. 

J i 


l'he  eddy  viscosity  v is  assumed  to  be  the  same  eddy  viscosity 
as  for  the  mean  flow.  Eliminating  r between  Equations  (5.8), 
(5.6)  and  (5.7),  we  find  that  the  coefficients  of  the  resulting 

^ du  d?u  . 
equations  are  functions  of  the  mean  flow  parameters  (u,— ,— — j ) 


and  v 


. In  order  to  solve  the  resulting  equations,  analytic 


expressions  of  mean  velocity  and  eddy  viscosity  have  to  be  spe- 
cified throughout  the  complete  turbulent  boundary  layer.  The  mean 
wind  profile  u(y)  was  measured  in  the  experiment,  but  the  phy- 
sical limiations  of  instruments  preclude  obtaining  experimental 
data  in  the  close  proximity  of  the  air-water  interface.  Con- 
sequently, the  analytic  expressions  of  mean  velocity  and  eddy 
viscosity  above  the  progressive  water  wave  have  to  be  assumed. 
Because  of  the  lack  of  experimental  data,  the  empirical  relation- 
ships describing  the  inner  and  o *.ir  regions  of  the  turbulent 
boundary  layer  over  a flat,  rigid  surface  are  used.  The  con- 
stants in  these  empirical  expressions  are  determined  with  the  best 
fit  to  the  experimentally  measured  mean  wind  profile.  Details  are 
presented  in  Sections  5.2.1  and  5.2.2 


5.2.1  Inner  Region 

Van  Driest 's  (1956)  eddy  viscosity  for  the  wall  region 
is  assumed  in  this  calculation  and  expressed  as  a non-dimensional 


. 


1 


F 


eddy  viscosity  based  on  v: 


1 + 4 K 


y+2£l  - exp(-y+/A+)j 


(5.9) 


where  y+  = yu*/v , K is  the  von  Karman  constant,  and  A+  is  the  para- 


meter characterizing  the  viscous  layer  thickness  on  the  y scale. 


A in  this  study  was  chosen  such  that  the  expressions  for  the  inner 


and  outer  layer  mean  velocity  profiles  will  overlap  at  y = 70. 
+ 


The  A values  for  different  wind  velocities  are  described  in 
Section  5.3.2. 

The  mean  velocity  in  the  wall  region  is  defined  by  using  Van 
Driest's  eddy  viscosity  expression  and  the  assumption  of  constant 
shear  stress  (Kline  et  al , 1968)  so  that 


u+(y+)  = 


f-T 

I 1 + 

J o 


2d  C 


1 + 4 k 2 7 2 1-  exp (- ?/A+) j ' 


0 1 y+  . 70 


(5.10) 


5.2.2  Outer  Region 

Coles'  (1956)  "low  of  the  wake"  expression  was  used 
for  the  outer  region  of  the  turbulent  boundary  layer.  The  velo- 
city correlation  in  the  outer  region  can  be  written  as 


U - u 


:i(x)  r 


-7n  ( y/6 ) 


+ — - 


1 + cos(iry/6) 


70v  / 1 


1.0 


(5.11) 
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where  6 is  the  boundary  layer  thickness,  70  v/(u*<$)  is  the  height 


corresponding  to  y+  = 70,  and  11  is  the  wake  parameter.  The  asymp- 
totic matching  condition  between  the  inner  and  outer  region  velo- 
city profiles  provides  (Stanford  Conference,  1968): 

211 

~ = —in  (<Su*/v)  + A + — (5.11a) 

* K 

U K 

* 

with  known  values  of  U , u , 6 and  A (obtained  by  least  square 

OO 

fitting  the  measured  mean  velocity  profile  to  u = — ■ ny  + A) 

Equation  5.11a  was  used  to  find  the  wake  parameter. 

The  eddy  viscosity  expression  was  chosen  such  that  it  is 
a constant  in  the  outer  region  and  approaches  <y  close  to  the  wall 
(Kline  et  al , 1968),  i.e., 

i,  2 du 

v = F 2(— ) — 

T V dy 

where  l is  the  mixing  length  determined,  from  experimental  results 
for  the  turbulent  boundary  layer  over  a flat  rigid  surface,  to  be 


> = 0.0858 6 tanh 


(_£_  ,V) 

''0.085  A; 


(5.12) 


and  F is  a correction  factor  determined  from  the  velocity  variation 
between  the  laminar  sub-layer  and  the  turbulent  region.  Using  the 
law  of  the  wall  proposed  by  Van  Driest  (1956) 

F = 1 - evr*  {-y  / A ) 
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A typical  comparison  between  the  measured  mean-velocity  pro- 
file and  the  analytical  expression  from  Equations  (5.10)  and  (5.11) 


is  shown  in  Figure  5.2.  The  analytic  expression  Equation  (5.11) 

representing  the  velocity  distribution  in  the  outer  region  repre-  j 

sents  the  experimental  data  quite  well.  (Note  that  the  lowest 

height  of  measurement  is  y+  = 200) . Although  the  experimental  data 

extends  smoothly  into  the  inner  region  represented  by  Equation 

(5.10),  the  adequacy  of  the  representation  (based  on  data  over  a 

rigid  surface)  cannot  be  confirmed  because  of  the  lack  of  detailed 

experimental  data  in  close  proximity  of  the  air-water  interface. 

Norris  and  Reynolds  (1975)  solved  the  perturbation  Equations 
(5.5  through  5.8)  in  the  transformed  coordinate  system.  Trans- 
formed Eddy  viscocity  Model  was  solved  numerically  by  EDVIS,  a 
FORTRAN  program  developed  by  Reynolds  and  Norris.  Details  relating 
to  transformation  of  the  governing  perturbation  equations  into  the 

~k  -k  -k 

(x  , y , t ) frame  and  the  computer  program  used  to  solve  the  trans- 
formed perturbation  equations  were  described  in  Norris  and  Reynolds 
(1975)  and  are  not  repeated  here.  In  the  transformed  coordinate 
system,  the  mean  velocity  and  eddy  viscocity  expressions  above  the 


wave  surface  were  assumed  to  be  the  same  as  the  expressions  as 

* 

described  in  Sections  5.2.1  and  5.2.2  with  y being  replaced  by  y . 
The  comparison  between  the  computed  results  and  the  experimental 
data  is  presented  in  section  5.6.  The  quasilaminar  model  compu- 
tation was  made  by  setting  the  eddy  viscosity  in  the  eddy  viscosity 
model  to  zero. 
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5.3  Turbulent  Kinetic  Energy  Model 

The  obvious  shortcoming  of  the  eddy  viscosity  model  lies  in 
the  implication  that  the  wave-induced  turbulent  Reynolds  stresses 
are  purely  a local  phenomenon  and  that  their  effect  at  each  point 
is  determined  by  the  character  of  the  mean  flow  at  that  point  alont 
The  flux  of  momentum  which  the  Reynolds  stresses  represent  is 
largely  the  result  of  mixing  between  regions  where  the  properties 
of  the  mean  flow  differ.  The  turbulence  which  accomplishes  the 
mixing  process  is  born,  convected,  diffused  and  dissipated  through 
viscous  action.  Consequently,  the  turbulent  structure  and  the 
mixing  depend  on  the  structure  of  the  flow  everywhere  and  are  not 
strictly  local  phenomena.  Prandtl  (1945)  recognized  this  diffi- 
culty and  suggested  that  the  turbulent  kinetic  energy  equation 
must  be  introduced. 

The  turbulent  kinetic  energy  equation  may  be  expressed  as 
(Norris  & Reynolds,  1975) 


9q'  _ 3q? 

+ u — = 2p-  2£)h 


3x . 3x . 
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(5.13) 


where  q'  = u.'u'.  Here, 


(P  - R.  . S.  . 

'•  tj  iJ 


(5.14) 


is  the  "turbulence  production". 


Is  the  "isotropic  dissipation"  (Reynolds,  1970). 
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Norris  and  Reynolds  (1975)  introduced  closure  assumptions  for 
(PJE).  the  pressure  velocity  term  and  the  triple  velocity  correlation 
term  as  follows: 

The  closure  model  for  the  turbulent  Reynolds  stresses  R..  is 

tj 

the  constitutive  equations 


R ..  = (—)  6 . . - 2 v S . . 
ij  3 iJ  T ij 


(5.15) 
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where  the  strain  rate  S..  = — ( — -1-  H ^) 

ij  2 ' . 
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In  Equation  (5.15),  the  turbulent  analogue  of  molecular  viscous 

stress  is  assumed  for  R,.  and  the  eddy  viscosity  v„,  = c e is 

ij  T i 

prescribed  as  a function  of  turbulent  kinetic  energy  and  a length 
scale  (Prandtl,  1945)  as  follows 


c q y/  v 

v = c e = c qt  (1  - e ) 

r 1 1 


(5.16) 


In  Equation  (5.16)  the  small  eddy  viscosity  characteristic  of 
the  viscous  sublayer  was  retained  by  the  near  wall  damping  factor 

c q y/  v 
3 

(1  - e ). 

In  order  to  minimize  the  difference  between  the  assumptions 
introduced  in  the  eddy  viscosity  and  turbulent  kinetic  energy 
models,  the  length  scale  l in  the  turbulent  kinetic  energy  model 


L&».  ' • \T  TgnrapMTmfr  ■ «sS£r  * •* 


was  chosen  to  have  the  same  characteristics  as  the  mixing  length 
specified  for  the  eddy  viscosity  model;  the  length  scale  ap- 
proaches > y in  the  wall  region  and  equals  a constant  (0.085)  in 
the  outer  region.  Hence, 


C = 0.0856  I 1 - (1 


- ,*>"] 


(5.17) 


where  n is  a constant  chosen  according  to  n x 0.085  = k , so  that 
in  the  wall  region,  the  length  scale  £ approaches  Ky.  The  pro- 


duction term|f  - R.  S for  the  special  case  of  a shearing  motion 
in  this  study  is  according  to  5.14 


<P 

dy 


(5.18) 
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The  dissipation  termjj  = v was  modeled  such  that  the 

Xj  j 

behavior  far  from  the  wall  is  controlled  by  the  large  eddy  expres- 
sionJJ»q!/£  (Kolmogorov,  1942)  and  near  the  wall  the  dissipation 
is  controlled  by  the  viscosity£)°c  Vq2/ii2.  The  combined  dissi- 
pations which  include  both  near  wall  and  far  from  the  wall  regions 
may  be  expressed  as 


£ = c -V  (i 


^ ^ * c£> 


(5.19) 


The  pressure  velocity  term  the  and  triple  velocity  correlation  term 
were  modeled  in  a gradient  diffusion  manner,  i.e.. 


(—  u.u.u.  + p u.  6.  )= (c  c— — ) 

2 l l 1 l i\  t 

>x.  J rlx.  ‘ <ix. 


(5.20) 


The  closure  assumptions  as  discussed  above  contain  five  free  con- 
stants, viz.,  c , c , c , c and  c . The  determination  of  these 

1/34 

constants  is  described  in  section  5.3.2. 

Norris  and  Reynolds  (1975)  adopted  the  turbulent  kinetic 
energy  equation  (5-13)  and  the  x-momentum  equation  together  with 
the  assumed  closure  condition  Equations  (5-16)  to  (5-20)  to  calcu- 
late the  mean  flow  quantities  for  a two-dimensional  turbulent  chan- 
nel flow.  Excellent  agreement  was  achieved  between  the  results  of 
the  mean  model  calculation  and  the  experimental  results  of  Hussian 
(1970)  and  I.aufer  (1951).  It  would  seem  natu  il  and  reasonable 
then  to  expect  that  these  closure  assumptions  may  be  applicable 
to  a two-dimensional  channel  flow  with  a wave-perturbed  motion. 
Consequently,  the  perturbed  turbulent  kinetic  energy  equation  was 
derived  by  formal  perturbation  of  Equation  (5-13)  with  assumed  gra- 
dient diffusion  term  and  with  Equation  (5-20)  for  the  pressure  velocity 
term  and  the  triple  velocity  correlation  term  in  Equation  (5-13). 

The  resulting  perturbation  turbulent  kinetic  energy  equation  is 


‘ ’ TV 4 


where  Q ‘ is  Che  mean  turbulent  kinetic  energy  for  perturbed  chan- 
nel flow  and  E is  the  eddy  viscoscity  function  for  the  mean  flow. 
When  the  formal  perturbation  technique  used  in  obtaining  Equa- 
tion (5.21)  is  applied  to  the  closure  condition  for  the  mean  flow, 
the  results  are  found  to  be 


where  Re  = OL/  v and  v =c  t.  Here,  in  contrast  to  the  eddy  viscosity 
1 ] 

model,  the  eddy  viscosity  is  allowed  to  have  perturbations  in  the 

turbulent  energy  model  as  shown  in  Equation  (5.23),  i.e.,  t.  = c q‘ 

6 

where  c is  a function  only  of  a mean  field  parameter.  Calcula- 

6 

tion  of  a perturbation  turbulent  energy  provides  a means  for  calcu- 
lating a wave  related  eddy  viscosity.  The  length  scale  L = 0.0 8 5 1 * 

- (1  - y*/6)njhas  been  assumed  as  a function  of  y"  and  the  func- 
tional relationship  is  the  same  as  Equation  (5.17)  with  y being 
replaced  by  y*.  Norris  and  Reynolds  solved  Equations  (5.5),  (5.6), 
(5.7),  and  (5.21)  together  with  the  closure  assumptions  (5.22) 
through  (5.25)  in  the  transformed  coordinate  system.  The  transformed 
Turbulent  Kinetic  Energy  Closure  Model  was  solved  numerically  by 
use  of  a FORTRAN  computer  program  developed  by  Reynolds  and  Norris. 
Again,  Che  mean  velocity  u(y*)  and  mean  eddy  viscosity  E(y  ) ex- 
pression above  the  wave  surface  in  our  study  were  assumed  to  be 
same  as  the  expressions  as  described  in  Sections  (5.2.1)  and  (5.2.2) 
with  y being  replaced  by  y*.  The  mean  turbulent  kinetic  energy 
expression  Q(y")  above  the  wave  surface  is  assumed  to  be  the  same 
as  that  of  a turbulent  boundary  layer  over  a flat  rigid  surface 

and  is  described  in  Section  (5.3.1).  The  five  free  constants  c , 

1 

c , c , c , and  c , in  Equations  (5.15),  (5.16),  (5.19)  and  (5.20) 

2 3 4 

are  dischussed  in  Section  (5.3.2). 

5 3.1  Mean  Turbulent  Kinetic  Energy 

As  shown  in  Equation  (5.21),  the  perturbation  turbulent 

energy  equation  has  non-constant  coefficients  determined  by  the 
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mean  turbulent  kinetic  energy  gradient.  Since  there  was  no  mean 

turbulent  energy  measurement  made  above  the  wavy  surface  in  this 

study,  the  analytic  expression  for  the  non-dimensional  mean  tur- 

■ "b  2 Q 2 

bulent  kinetic  energy  (based  on  u*)  Q = —yrr  was  obtained  by  a 

u 

curve  fit  of  Klebanoff's  (1954)  turbulent  boundary  layer  experi- 
mental results  (above  a solid  rigid  surface)  with  the  expression 

Q+‘  = By+(2_d)  exp  - d2(tny+  - b)  ' + F (5.26) 

where  B = exp  (-0.012),  d = 1.115,  b = 2.94  and  F = 0.0001.  Figure 
(5.3)  shows  the  comparison  between  the  experimental  data  and  the 
curve  fit  of  Equation  (5.2b). 


5.3.2  The  Model  Constan ts 

From  Equations  5.15,  5.16,  5.19  and  5.20,  the  turbulent 

kinetic  energy  model  contains  five  free  constants  c , c , c , c 

12  3 4 

and  c which  need  to  be  determined.  From  several  known  characteris- 

5 

tics  of  turbulent  boundary  layer  flows  these  constants  can  be 
either  closely  or  precisely  determined.  First,  in  the  logarithmic 
section  of  the  layer  the  slope  of  the  mean  velocitv  profile  in  semi- 
logarithmic  wall  coordinates  is  1/-  . production  and  dissipation 
are  approximately  equal  and  the  u ssipation  at  the  wall  is  a non- 
zero finite  value.  The  constants  chosen  should  reflect  these  facts 

+ 

and  experimentally  measured  A values  and  q'  profiles.  Details 
of  the  derivation  of  the  constants  are  given  in  Norris  and  Reynolds 
(1975)  and  Achary  and  Reynolds  (1975)  and  are  not  repeated  here. 


- 


The  main  features  in  obtaining  those  constants  and  the  constants 


used  in  this  study  are  as  follows: 

L.  Norris  and  Reynolds  integrated  the  momentum  equation  for 

the  turbulent  channel  flow  in  wall  coordinates  and  obtained 

c =1  /q 
1 

by  using  the  slope  of  1/k  as  noted  above.  Based  on  Laufer's 
(1951)  channel  flow  data, 

5.0  < q+  ' < 8.0 

in  the  logarithmic  region  extending  from  30  < y+  < 300.  They 
+ 2 

chose  a value  of  q = 6.5  (c  close  to  k is  a good  number  to  choose) 

1 

■f"  2 

and  obtained  c = 0.39.  In  this  study  we  chose  q = 5.95  and 

1 

obtained  c = 0.41. 

1 

2.  Norris  and  Reynolds  assumed  that  the  turbulent  diffusivity 

for  turbulent  kinetic  energy  and  momentum  should  be  approximately 

equal  to  each  other,  therefore  c = c . Their  calculations  showed 

2 1 

that  the  model  is  not  sensitive  to  the  c value  chosen  (they  varied 

2 

the  c value  from  0.5  c|  to  2 c^. 

3.  Considering  that  the  eddy  viscosity  damping  factor 

+ + 

-c  q y 

(1  - e f ) in  Equation  (5.16)  and  the  Van  Driest's  damping 

factor  (1  - e ^ ^ ) play  the  same  role,  clearly  c = l/q+A+. 

3 
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This  constant  reflects  experimentally  measured  A and  q values. 

As  described  in  Section  5.2.1(the  A+  in  this  study  was  obtained 
by  matching  the  experimentally  obtained  mean  velocity  profile  in 
the  outer  region  (curve  fit  with  the  Equation  5.11)  with  the  as- 
sumed inner  region  velocity  expression  of  Equation  5.10.  Kays  (1971) 


In  this 


showed  that  A = 26  for  the  flat  plat  boundary  layer. 

study,  when  U /C  = 0.883,  1.115  and  2.283,  the  corresponding  A 

values  are  8.42,  7.55  and  6.18,  respectively.  This  suggests 

that  the  velocity  profile  above  the  progressive  water  wave  has 

a thinner  sublayer  as  compared  with  that  above  the  flat  plate. 

In  determining  the  q+  value,  we  used  Klebanoff's  experimental 

data  for  turbulent  boundary  layer  flow  that  in  the  region  (20  _ 
y+  _ 40)  strongly  affected  by  the  eddy  viscoscity  damping  term, 

q+  = 12,  therefore  q+  = 3.464.  In  Norris  and  Reynold's  study, 

they  used  Lanfer's  channel  flow  data  and  chose  q+  = 8. 

4.  From  experimental  data,  it  is  known  that  in  the  log  region 

production  and  dissipation  approximately  balance  each  other.  If 

production  and  dissipation  are  equated 

3 

c = c 

4 1 


5.  Norris  and  Reynolds  used  the  facts  that  at  the  wall  the 

i 

dissipation  is  a finite  non-zero  value  and  turbulent  production 


and  eddy  viscosity  are  zero  to  determine  from  Equation  (5.19)  and 


the  turbulent  energy  equation  that 

I 2. 

c = k /c 

i ■ 6 4 

Therefore,  a typical  set  of  values  of  the  model  constants 


based  on  the  above  derivations  is 


8.42  when 


c =i/^+a+)  = 0.03428  (q+  = 3.464  and  A+  = 

3 

U 

= 0.883) 

c = c 3 = 0.06892 

4 1 

c = k^/c  = 2.4  39 

5 4 

5 . 4 Boundary  Conditions 
The  boundary  conditions  for  a channel  flow  with  progressive 

waves  at  the  lower  channel  boundary  are  as  follows.  Because 

all  fluid  velocities  at  the  flat  stationary  upper  wall  are  zero, 

~ - ~ 2 

the  no  slip  condition  requires,  u = v = q =0  there.  The 
lower  boundary  is  a progressive  wave  and  the  kinematic 
boundary  condition  requires  that  u = CKn  and  v = -iKCn  at  the  wave 
surface.  Because  of  the  effect  of  the  wind,  small  ripples  are 
generated  on  the  surface  of  the  mechanically-generated  waves. 
Therefore,  the  turbulence  energy  need  not  vanish  at  the  wave  sur- 
face. In  considering  the  undertainties  regarding  the  wind  field 
and  wave  field  interaction  and  the  sensitivity  of  the  wave- 
induced  turbulent  energy  boundary  condition  to  the  computation 
(Long,  1970)  , we  assumed  that  all  background  turbulence  velo- 
city components  are  zero  at  the  air-water  interface,  i.e., 
q = f'  at  wave  surface. 

5 . 5 Model  Solutions  and  Compari sons  with  the  Data 

Norris  and  Reynolds  solved  the  perturbation  equations  in  the 
transformed  coordinate  system  described  above  so  that  the  boundary 
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conditions  could  be  applied  directly  at  the  wave  surface.  Since 
the  wave  follower  measurement  is  neither  in  the  (x,  y,  t)  nor 

^ JL 

in  the  t ransformed (x  , y , t ) coordinate  systems,  the  fixed 
probe  measurements  were  compared  with  the  model  predictions.  After 
the  computed  results  of  wave-induced  velocities  and  turbulent 
Reynolds  stresses  were  obtained  in  the  transformed  coordinate 
system,  the  following  steps  were  used  to  obtain  the  corresponding 
results  from  the  fixed  probe  measurements. 

(a)  Obtain  the  values  of  |u  (y*)|  ,|v(y*)  and  r _(y*) 
from  the  modified  Norris  and  Reynolds  (1975)computer  program  as 


a functions  of  y“  and  phase  angle  (a  total  of  11  phase  angles  with 
30°  increment)  with  respect  to  wave.  Based  on  the  coordinate  trans- 
formation their  corresponding  values,  i.e.  |u(y)|,  |v(y)|  and 

Jr . j (y)J  were  calculated.  Note  that  |u(y)|  is  related  to  |u(y*)| 
by 

|u(y)|  = |u(y)|  - ^ * nf  + 0 (n2) 
where  f is  given  by  Equation  (5.4) 

(b)  Pick  up  the  results  of  | u(y)  | , | v(y)|  and  | r„(y)|  from 
(a)  for  a particular  value  of  y over  different  phases  of  the  wave. 

(c)  Lease-square  curve  fit  the  results  of  (b)  with  a sinu- 


soidal curve  to  obtain  the  magnitude  and  phase  of  u. 


v and  r . . 

ij 


for  a fixed  height  above  the  mean  water  surface. 


(d)  Repeat  (b)  and  (c)  to  obtain  the  magnitude  and  phase 


moot ) . The  model  computations  were  carried  out  for  three  different 


mean  flov  conditions,  and  the  results  are  presented  below. 


5.5.1  No  Critical  Layer  Case  (U  /C  = 0.883) 


Li 


Saeger  and  Reynolds  (1971)  found  that  their  turbulent 
models  were  not  adequate  for  the  flowfield  with  a critical  layer 
and  better  agreements  were  achieved  for  upstream  running  waves,  i.e., 
cases  with  no  critical  layer  in  the  flowfield.  Therefore,  it 
is  interesting  to  see  whether  similar  conclusions  can  be  reached 
here. 

In  Figure  4.3  the  trends  of  computed  |u  1 and  |v  j are  similar 
to  those  of  the  experimental  data  for  both  the  eddy  viscosity 
and  turbulent  energy  models.  However,  the  computed  values  of  |u  j 
and  |v  | are  consistently  higher  than  the  experimental  data.  In 
Figure  (4.4),  the  computed  and  experimental  phase  lags  agree  well. 

The  comparisons  of  the  magnitude  of  r _ are  shown  in  Figures  4.13a, 
b,  c.  Both  computer  models  predicted  the  same  order  of  magnitude 
for  all  three  components  or  r„,  while  the  experimental  results 
show  order  of  magnitude  difference  among  the  three  components  of 
r in  the  region  close  to  the  wave  surface.  The  comparison  of 
the  phase  lag  of  r is  shown  in  Figure  4.14a,  b,  c.  Clearly,  both 
models  fail  to  predict  both  the  magnitude  and  the  phase  angle 
variation  as  a function  of  non-dimensional  distance  Kh. 
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5.5.2 


Thick  Critical  Layer  Case  (U  /C  = 1.115) 


The  critical  layer  height  is  12.38  cm  above  the  mean 
water  surface.  As  shown  in  Figure  4.6,  the  experimental  data  show 
a sharp  phase  shift  of  the  u-component  near  the  critical  layer 
height.  Such  a behavior  may  be  used  as  an  indication  of  the  capa- 
bility of  prediction  of  various  closure  models.  The  quasilaminar 
model  was  included  also  in  this  case  by  setting  the  eddy  viscosity 
in  the  eddy  viscosity  model  to  zero  to  reveal  the  sensitivity  of 
the  r terms.  Figure  4.5  shows  the  comparison  for  the  magnitude 
of  wave-induced  velocities.  Again,  all  the  computational  results 
are  larger  than  those  of  the  experimental  data.  However,  the 
difference  between  models  is  small.  Figure  4.6  is  the  comparison 
for  the  phase  lags  of  wave-induced  velocities.  The  computed  results 
agreed  well  with  the  experimental  data  close  to  the  wave  surface, 
but  fail  to  follow  the  trend  of  the  experimental  data  away  from 
the  surface.  The  eddy  viscosity  and  quasilaminar  models  yield  the 
same  trend  for  the  phase  angle  of  the  v-component,  while  the  eddy 
viscosity  and  turbulent  energy  models  yield  the  same  trend  for  the 
phase  angle  of  u-component.  None  of  the  models  predict  the  large 
phase  variation  of  u near  the  critical  layer.  The  comparisons  of 
r between  experimental  data  and  computed  results  are  shown  in 
Figures  4.15a,  b,  c.  Similar  to  the  non-critical  layer  case,  the 
models  yield  the  same  order  of  magnitude  on  all  |r^.  (while  the 
experimental  data  indicate  differences  in  their  magnitude,  !r 


being  the  largest  compared  to  both  r 1 and  'r  . The  corres- 
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ponding  phase  lags  of  r terms  are  shown  in  Figure  4.16a,  b,  c. 

The  phase  lags  of  r predicted  by  the  turbulent  energy  model  agree 
1 1 

well  with  the  experimental  data,  while  the  eddy  viscosity  model  pre- 
dicted higher  values.  Both  models  fail  to  predict  the  phase  lags 

of  r and  r variation  as  function  of  Kh . 
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5.5.3  Thin  Critical  Layer  Case  (U  /C  = 2.283) 


In  this  case,  the  Reynolds  number  is  the  largest  among 
all  the  experimental  runs  conducted.  Because  of  the  high  wind 
speed  the  critical  layer  was  inaccessible  to  the  probe.  Only  the 
eddy  viscosity  model  was  analyzed  in  this  case.  Figure  4.11  shows 
the  comparison  between  the  experimental  and  computed  magnitude 
of  wave-induced  velocities.  Although  the  trends  agree  well  with 
the  experimental  data,  the  computed  results  are  higher  than 
those  of  the  experiments.  Figure  4.12  shows  the  phase  lag  com- 
parison for  the  wave-induced  velocities  and  the  trends  compare 
favorably  with  those  of  the  experiment.  As  shown  in  Figure  4.21a, 
b,  c and  4.22a,  b,  c,  the  agreement  between  the  computational  and 
experimental  results  are  poor  for  both  the  magnitude  and  phase  of 
wave- induced  turbulent  Reynolds  stresses.  The  model  cannot  pre- 
dict the  order  of  magnitude  difference  among  r..  in  the  region 
close  to  the  interface  and  cannot  predict  the  trend  of  the  phase 
lags  of  r_  particularly  when  far  from  the  interface. 


A 


5 . 6 Discussion 

In  Norris  and  Reynolds  (1975)  the  ratio  of  wave  length  to 
wave  amplitude  is  281,  the  Reynolds  number  is  13,958,  and  C/l  for 

the  downstream  running  wave  is  in  the  range  0 1 C/U ^ - 0.4.  In 
this  study,  the  wave  length  to  amplitude  ratio  is  62,  the  Reynolds 
number  varies  from  43,670  to  112,862  and  C/U^  is  in  the  range 
0.43  _ C/U^  2 1.13.  Norris  and  Reynolds  (1975)  (perturbation  pres- 
sure calculation  in  channel  flow  of  moderate  Reynolds  number)  indi- 
cates that  in  the  case  of  a long  wave,  the  use  of  coordinate  trans- 
formation to  apply  the  boundary  conditions  directly  at  the  wave 
surface  seems  crucial  and  the  choice  of  Reynolds  stress  model  is 
irrelevant.  In  our  study  the  computed  results  from  the  models  of 
Norris  and  Reynolds  do  not  agree  with  the  measured  wave-induced 
perturbations.  The  results  are  conclusive,  but  the  reader  should 
note  the  following: 

(a)  The  coefficients  of  the  wave-perturbed  momentum  and 
energy  equations  are  functions  of  the  mean  velocity  and  turbulent 
energy  across  the  turbulent  boundary  layer.  Due  to  the  fact  that 
the  experimental  data  in  close  proximity  to  the  air  water  inter- 
face could  not  be  obtained,  the  mean  velocity  profile  (below  y+  = 
200)  and  the  turbulent  energy  above  the  water  waves  used  to  cal- 
culate the  coefficients  were  adapted  from  the  data  of  turbulent 
boundary  layer  flow  over  a flat  rigid  surface.  Also,  the  require- 
ment that  the  velocity  be  continuous  across  the  air-water  inter- 
face restricts  the  mean  velocity  at  the  mobile  wave  surface  to 
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be  equal  to  the  mean  current.  With  the  adapted  mean  velocity  pro- 
file, the  velocity  is  zero  at  the  air-water  interface.  Davis 
pointed  out  in  both  of  his  studies  (1970,  1972)  that  the  model 
calculation  is  very  sensitive  to  the  mean  velocity  profile  close 
to  the  air-water  interface. 

(b)  The  uncertainties  inherent  in  the  choice  of  the  proper 
boundary  condition  to  apply  at  the  wave  surface  are  the  following: 
The  water  surface  is  roughened  by  the  turbulent  shear  flow,  hence 
the  turbulent  energy  need  not  vanish  at  the  interface  and  it  is 
possible  that  there  exists  an  oscillating  component  of  turbulent 
energy  which  is  correlated  with  the  wave.  Moreover,  the  eddy 
viscosity  which  used  to  represent  the  turbulent  Reynolds  stress 
at  the  interface  need  not  vanish  at  the  roughened  wave  surface. 

In  this  study,  both  the  eddy  viscosity  and  turbulent  energy  are 
assumed  to  be  zero  at  the  wave  surface.  Long  (1971)  used  the  non- 
zero boundary  condition  of  turbulent  energy  at  the  wave  surface 
and  carried  out  several  calculations;  the  computed  induced  pres- 
sure amplitude  was  too  large  by  several  orders  of  magnitude.  His 
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large  as  the  mean  turbulent  kinetic  energies.  Therefore,  the 
linear  assumption  used  in  their  and  in  our  studies  for  the 
model  computation  might  not  be  applicatble. 


6.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  overall  objectives  of  this  investigation  wen-  largely 
accomplished.  A comprehensive  and  reliable  set  it  experimental 
data  on  both  the  magnitude  and  the  phase  of  wa'  e-induced  veloci- 
ties and  turbulent  Reynolds  stress  was  acquired  across  a turbu- 
lent boundary  layer  above  a small  amplitude  1 Hz  progressive 
water  wave.  The  data  reveals  the  importance  of  the  wave-induced 
turbulent  Reynolds  stress  in  the  momentum  transfer  from  wind 
to  water  waves.  Tn  addition,  further  numerical  calculations  based 
on  various  closure  models  for  predicting  the  wave-induced  pertur- 
bations in  a turbulent  shear  flow  should  be  postponed  until  more 
basic  experimental  data  on  the  structure  of  turbulence  above  the 
air-water  interface  is  available. 

The  specific  conclusions  derived  from  this  study  and  the 
recommendations  made  for  the  future  study  are  as  follows: 

(1)  The  quasilaminar  model  (r  = 0)  is  not  adequate  for 
description  of  water-wave  generation  by  the  wind.  This  study 
indicates  that  a significant  contribution  to  the  total  momen- 
tum transfer  from  wind  to  wave  is  derived  from  the  r. . terms 

ij 

and  from  the  region  close  to  the  air-water  interface. 

(2)  It  is  clear  from  the  discussion  of  Section  5.6  that 
the  experimental  data  for  mean  velocity,  turbulent  energy,  and 
turbulent  Reynolds  stresses  above  the  mobile  deformable  wave 
surface  are  important  to  clarify  the  adequacy  of  the  eddy  vis- 
cosity and  turbulent  energy  models  used  in  this  study.  Both  the 
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fixed  and  wave  following  frame  measurements  as  used  in  this  study 
cannot  obtain  the  data  at  the  air-water  interface.  The  develop- 
ment of  the  Laser  Dopier  Anemometer  Measurement  Technique  may 
allow  us  to  obtain  data  in  close  proximity  to  the  air-water  inter- 
face . 

(3)  The  effects  of  drift  current  and  the  superposed  capillary 
waves,  which  ride  on  the  progressive  water  waves,  on  the  turbulent 
structure  should  be  studied  together  with  the  detailed  wave 
height  signal  to  reveal  better  information  for  determining  the 
empirical  function  for  the  closure  models.  Uzkard  and  Reynolds 
(19b7)  studied  the  effect  of  a moving  rigid  wall  on  the  structure 
of  a turbulent  boundary  layer.  As  indicated  in  their  experimental 
results,  the  turbulent  fluctuations  were  suppressed.  If  the 
current  at  the  air-water  interface  is  comparable  to  the  movement 

of  a solid  boundary,  the  turbulent  structure  above  the  air-water 
interface  would  be  different  from  that  above  the  flat  rigid  sur- 
face. Also,  the  turbulent  energy  and  Reynolds  stresses  need  not 
vanish  at  the  interface  because  the  water  surface  is  roughened  by 
the  turbulent  shear  flow. 

(4)  The  magnitude  of  the  mean  quantities,  e.g.,  velocity,  eddy 
viscosity  and  turbulent  kinetic  energy,  are  assumed  to  be  con- 

k 

stant  along  a constant  y coordinate  in  the  model  computation  of 
this  study.  In  other  words,  the  functional  relationships  in  the 
transformed  coordinate  system  are  assumed  to  be  the  same  as  those 
in  the  unt ransformed  coordinate  system  except  for  the  coordinate 

k 

y being  replaced  by  y . This  assumption  seems  to  be  physically 
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realistic  because  values  adjust  themselves  in  accordance  with  the 
passing  waves  instead  of  being  constant  in  the  untransformed  co- 
ordinate system.  The  above  described  assumptions  should  be  veri- 
fied experimentally. 

(5)  According  to  the  study  of  Acharya  and  Reynolds  (1975)  , 
the  turbulent  energy  model  cannot  predict  the  transfer  of  energy 
between  the  different  components  of  r„  and  adequate  modeling  of 
the  pressure  strain  terms  in  the  dynamic  equations  of  r^  is  cru- 
cial if  success  is  to  be  achieved.  Therefore,  modeling  the  flow 
by  using  the  turbulent  Reynolds  stresses  equations  could  be  of 
great  interest  and  utility. 
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APPENDIX  I 


UNCERTAINTY  ANALYSIS 


When  a result,  R,  Is  computed  from  several  different  inputs,  each 
of  which  has  an  uncertainty,  then  the  value  of  the  result  is  uncertain. 
If  the  uncertainties  in  the  inputs  are  bounded  then  we  can  describe 
the  greatest  possible  excursion  in  R which  could  occur  due  to  the 
least  favorable  combination  of  the  largest  excursions  in  each  input. 

Kline  and  McClintock  (1953)  have  shown  that  the  way  to  combine 
the  individual  effect  which  most  nearly  preserves  the  true  statistical 
probability  is  to  use  a root-sum-square  addition: 
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Where  the  <5x^'s  are  the  independent  random  uncertainties  with  normal 
distributions  upon  which  R depends  and  6R  is  the  uncertainty  interval 
in  R,  expressed  at  the  same  confidence  level  as  was  used  in  describing 
the  inputs.  Note  that  only  the  largest  error  and  terms  of  like  order 
of  magnitude  will  significantly  contribute  to  the  uncertainty  in  R. 

The  following  uncertainties  were  used  in  the  uncertainty  calcu- 
lation. 


Temperature 
Pressure  difference 
Calibration  of  hot  film 
Calibration  of  wave  height  gauge 


±0.006  cm  of  0.82  S.  G.  fluid 
±4%  of  calibration  curve 
±1%  of  mean 
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The  resulting  maximum  uncertainties  of  wave  induced  perturbations 


computed  by  data  reduction  program  and  Equation  (A-l)  are 


Magnitude  of 
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+ 4% 
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Phase  of  r . . 
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0.  . 

Besides  the  single-sample  error  estimate  as  described  above,  the 
uncertainty  in  the  spectral  result  is  estimated  by 


where  N is  the  number  of  spectra  used  for  the  average,  and  £ is  the 
error  in  the  spectra.  The  resulting  uncertainty  of  wave  induced  per- 
turbation for  60  spectra  average  is 


Magnitude  of  u.  and  r..  ± 13%  of  u and  r 

1 ij  1 i‘  1 i j I 

Phase  of  u.  and  r..  ± 8°  of  0 and  9 

i il  i ij 

The  statistical  uncertainty  due  to  the  use  of  a finite  data  record 
is  much  greater  than  the  uncertainty  due  to  instrument  errors. 
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APPENDIX  II 


EXPERIMENTAL  DATA  1 


A.  Fixed 
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K - ^ 

= 4 .019  rad/m 

for  all  cases. 
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Figure  2.1  Cross  Correlation  and  Decomposition 
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Figure  3.2  Schematic  of  Probe  Carriage 
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Wave  Follower  Installation  Outside  the 
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Figure  3.5  Pressure  Transducer  Calibration  Curves 
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Figure  3.6  Capacitance  Type  Wave  Gauge 
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Figure  3.10  Typical  Coherence  Spectrum  of 
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Figure  4.6  Phase  Distributions  of  u and  v and 
Computed  Results  for  U^/C  = 1.115 
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Figure  4.9  Magnitude  Distributions  of  u and  v 

for  U /C  - 1.836 
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Figure  4.11  Magnitude  Distributions  of  u and  v and 
Computed  Results  for  U^/C  = 2.283 


Figure  4.12  Phase  Distributions  of  u and  v and 
Computed  Results  for  U^/C  = 2.283 
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Figure  4.13a  Magnitude  Distribution  of  r and 

1 1 

Computed  Results  for  iWC  = 0.883 
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Figure  4.13b  Magnitude  Distribution  of  r^  and 
Computed  Results  for  U^/C  = 0.883 
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Figure  4.14c  Phase  Distribution  of  r and  Computed 

Results  for  U /C  = 0.883 
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Figure  4.15a  Magnitude  Distribution  of  r 
Computed  Results  for  U^/C  = 


(aajfiap)  to-)  noMd 


Figure  4.16c  Phase  Distribution  of  r and  Computed 

1 2 

Results  for  U /C  = 1.115 
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Figure  4.18  Phase  Distributions  of  r 
for  U /C  = 1.34 
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Figure  4.27  Magnitude  Distribution  of  r : 

1 2 

Wave  Follower  Results 
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Figure  4.28a 
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Wave  Follower  Results 
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Figure  4.28b  Phase  Distributions  of  r^  : 
Wave  Follower  Results 
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Figure  5.1  Lines  of  Constant 


Figure  5.3  Comparison  of  Turbulent  Kinetic  Energy 


